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ABSTRACT 


Weightlessness  is  the  major  eontributing  faetor  behind  the  degradation  of  bone 
mass,  musele  tone,  and  aerobie  eapaeity  during  long-term  spaee  missions.  With  the  loss 
of  bone  mass  progressing  at  up  to  two  pereent  per  month,  long  duration  and 
interplanetary  missions  shall  remain  the  sole  duty  of  robotie  explorers  until  suffieient 
eountermeasures  are  developed. 

Several  eountermeasures  are  either  in  use,  or  under  development  to  alleviate  this 
problem.  Exereise  is  eurrently  used  to  reduee  the  severity  of  bone  loss  and  musele 
atrophy.  Exereise  has  proven  ineffeetive  despite  the  faet  two  hours  of  daily  exereise 
together  with  elaborate  apparatus  have  been  devoted  to  simulating  the  load  of  Earth’s 
gravity.  Drug  therapy  and  other,  more  exotie,  eountermeasures  are  also  under 
eonsideration,  but  the  side-effeets  of  these  other  treatments  and  the  faet  that  they  do  not 
direetly  address  the  root  eause  of  the  negative  effeets  of  weightlessness  means  that  they 
may  only  reduee,  not  eure,  those  problems.  Only  artifieial  gravity  addresses  the  root 
eause,  weightlessness  itself 

This  thesis  addresses  the  need  to  balanee  the  effeets  of  Coriolis  on  human 
disorientation  with  the  engineering  eosts  of  eonstrueting  a  eentrifuge  for  human 
oeeupation  in  spaee. 
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I.  INTRODUCTION 


A,  NEGATIVE  EFFECTS  OF  LONG  TERM  SPACE  FLIGHT 

Long  Term  Space  Flight  carries  with  it  numerous  negative  effects  upon  the  human 
body.  The  micro-gravity  environment  experienced  by  astronauts  in  space  imparts  a 
weakening  of  muscles,  bone  structure,  and  a  redistribution  of  bodily  fluids.  Effective 
countermeasures  have  been  developed  against  fluid  redistribution  and  muscular  atrophy 
to  mitigate  their  adverse  effects  upon  the  astronaut’s  return  to  Earth.  However,  effective 
countermeasures  have  not  been  developed  to  overcome  the  weakening  of  the  body’s 
bones  in  space  as  shown  in  Eigure  l.i 


Normal  vertebrae  Vertebrae  suffering 
from  osteoporosis 

Figure  1:  A  Comparison  of  Normal  Vertebrae  and  Vertebrae  Effected  By 

Osteoporosis^ 

1,  Bone  Loss  Associated  with  Space  Flight 

As  a  result,  bone  loss  is  one  of  the  greatest  obstacles  to  Long  Duration  Manned 
Space  Flight.  Osteoporosis,  a  well-publicized  disease  that  is  normally  associated  with 


1  Julie  Moberly,  “Human  Physiology  Research  and  the  ISS:  Staying  Fit  Along  the  Journey,”  Space 
Research.  March  2002,  pp.  6-11  and  25 

2  Scott  M.  Smith,  “Calcium  Kinetics  During  Spaceflight,”  19  April  2002,  Nutritional  Biochemistry 
Laboratory.  NASA  Johnson  Space  Center, 

http://www.isc.nasa.gov/sa/sd/facilitv/labs/nutritionalbl/currentproi/cal  kin.htm 
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advanced  age  and  menopause,  oecurs  when  the  body’s  natural  process  of  bone  renewal, 
ealled  remodeling,  breaks  down.  Remodeling  is  effected  by  a  eomplex  interplay  of  many 
hormones,  the  pull  of  museles,  and  gravity^.  Although  there  are  no  documented  eases  of 
Osteoporosis  caused  by  space  flight,  bone  loss  occurs  at  a  significant  rate  in  spaee: 

"When  you  remove  gravitational  loading,  bones  no  longer  sense  the 
stresses  and  strains  that  are  normally  experieneed  here  on  Earth.  As  a 
result,  astronauts  are  subjeeted  to  an  aeeelerated  rate  of  bone  loss,  losing 
between  a  half  of  1  pereent  and  2  pereent  of  their  bone  mass  per  month,"^. 

“Weakening  of  the  bones  due  to  the  progressive  loss  of  bone  mass  is  a 
potentially  serious  side-effect  of  extended  spaeeflight.  Studies  of 
cosmonauts  and  astronauts  who  spent  many  months  on  spaee  station  Mir 
revealed  that  spaee  travelers  can  lose  (on  average)  1  to  2  pereent  of  bone 
mass  eaeh  month. 

Without  some  means  of  preventing  or  redueing  space  flight  bone  loss,  astronauts 
on  a  theoretieal  two  year  mission  to  Mars  eould  lose  up  to  24%  of  their  bone  mass. 

2,  Methods  to  Counter  Bone  Loss  During  Space  Flight 

A  complex  interplay  of  physiology,  exercise,  and  environment  serve  to  maintain 
the  body’s  fluid  distribution,  muscle  tone,  and  bone  mass.  Because  our  bodies  evolved 
on  Earth,  they  are  designed  to  resist  the  eonstant  pull  of  gravity.  In  spaee  our  body’s 
natural  process  of  maintenance  and  repair  works  to  expel  fluids  displaeed  by  the  absence 
of  gravity,  muscles  atrophy  from  disuse,  and  the  remodeling  proeess  slowly  degrades 
bone  mass  due  to  the  absence  of  gravity®.  There  are  three  basie  means  of  eounteraeting 
the  detrimental  effeets  of  spaee  flight:  exercise,  drug  therapy,  and  artifieial  gravity. 


3  “Boning  Up  on  Osteoporosis”,  September  1996  (Revised  September  2001),  FDA  Consumer.  U.S. 
Food  and  Drug  Administration,  25  June  2002,  http://www.fda.gov/fdae/features/796  bone.html 

4  Ted  Bateman,  BioServe  Spaee  Teehnologies,  Prineipal  Investigator  and  Direetor  of  Biomedieal 
Researeh,  http://spaeereseareh.nasa.gov/general  info/issphvsiology.html 

®  Doug  Hullander,  Patriek  L.  Barry,  “Spaee  Bones,”  1  Oet.  2001,  Seienee@NASA.  NASA  Marshall 
Spaee  Flight  Center,  14  Sept.  2002,  http  ://seienee.nasa.  gov/headhnes/v200 1/astO  1  oet  1  .htm 

®  Moberly 
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a.  Exercise 

Aerobic  exercise  can  reduce  the  effect  of  some  of  the  physiological 
problems  associated  with  space  flight.  Astronauts  receive  individualized  exercise 
prescriptions  before  and  during  their  missions  that  are  designed  to  maintain  their  aerobic 
capacity,  bone  density,  and  muscle  mass  as  much  as  possible^.  However,  reducing  the 
effects  of  muscle  atrophy  currently  requires  about  two  hours  of  exercise  per  day  with  the 
aid  of  ‘exotic  devices’  to  reproduce  the  effects  of  Earth’s  gravity.  Unfortunately, 
exercise  has  proven  ineffective  as  a  countermeasure  to  bone  and  muscle  loss* *. 


Figure  2:  An  Astronaut  Exercising  on  a  Shuttle  Treadmill, 
Image  courtesy  of  NASA  Johnson  Space  Center^ 


b.  Drug  Therapy 

The  remodeling  process  manages  the  development  and  maintenance  of 
bone  health.  In  mature  adults,  the  remodeling  process  maintains  the  structure  of  normal 
bone  tissue  (pictured  in  Figure  3)  through  an  equilibrium  between  the  activity  of 
osteoblasts  (bone  forming  cells)  and  osteoclasts  (bone  absorbing  cells).  In  the  astronauts 


^  Moberly,  pp.  6  and  8 

*  Patrick  L.  Barry,  “Good  Vibrations,”  Science@,NASA.  2  Nov.  2001,  NASA  Marshall  Space  Flight 
Center,  14  Sept.  2002,  http://science.nasa.gov/headlines/v2001/ast02nov  l.htm 

^  Exercise  In  Space.  22  April  1998,  Exercise  Countermeasures  Project,  NASA  Johnson  Space  Center, 
14  Sept.  2002,  http://www.isc.nasa.gov/sa/sd/sd3/exl/spacephoto.htm 
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and  the  elderly,  osteoblast  activity  is  inhibited  while  osteoclast  activity  remains  constant. 
The  resulting  loss  in  bone  mass  eventually  leads  to  osteoporosis. 

Drug  therapy  could  potentially  reduce  the  amount  of  bone  loss 
experienced  by  astronauts  during  space  flight.  Alendronate  Sodium  (Fosamax),  for 
example,  inhibits  the  resorption  of  bone  by  inhibiting  osteoclast  activity,  reducing  or 
potentially  reversing  the  progression  of  osteoporosis.  However,  this  drug  is  poorly 
absorbed  following  oral  administration.  It  must  be  administered  separately  from  other 
medications  and  cannot  be  administered  with  food  or  caffeine.  It  also  increases  the  risk 
of  GI  problems  if  any  Non-Steroidal  Anti-inflammatory  Drugs  (NSAIDs)  are  usedif 


Figure  3:  The  Structure  Of  Normal  Bone  Tissue,  Image  Courtesy  of  NASA  Quest^^ 

In  addition  to  Alendronate  Sodium,  there  are  several  other  drugs  on  the 
market  that  are  used  to  treat  Osteporosis.  Risedronate  (Actonel),  like  Alendronate 
Sodium,  alters  the  remodeling  process  by  inhibiting  osteoclast  activity,  with  a  similar  list 
of  precautions  and  side  effects. 1 3  Raloxifene  (Evista)  works  like  estrogen  in  post- 

10  Moberly,  p.  10 

1 1  Judith  H.  Deglin,  April  H.  Vallerand,  Davis’s  Drug  Guide  for  Nurses.  5*  Ed.,  Philadelphia:  F.  A. 
Davis,  1997,  pp.  21-23 

13  Barry 

13  “Risendronate,”  WebMD.  2  June  2000,  Ver.  4.01,  Multum  Information  Services,  14  Sept.  2002, 
http://mv.webmd.eom/content/article/4046.1450 
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menopausal  women  to  prevent  osteoporosis. Finally,  Caleitonin  (Caleimar,  Miaealein) 
is  a  naturally  oeeurring  hormone  that  helps  to  regulate  calcium  levels  and  is  involved  in 
the  bone  building  part  of  the  remodeling  process.!^ 

Of  the  currently  marketed  drugs,  only  calcitonin  would  be  useful  to 
astronauts.  Raloxifene  is  intended  for  women  who  choose  not  to  take  estrogen  or  other 
medications.!®  Alendornate  Sodium  and  Risedronate  is  intended  for  patients  that  can 
remain  upright  for  at  least  30  minutes  and  then  eat  after  taking  the  medication  to  prevent 
heartbum!^’!*’!^,  (impossible  in  a  micro-gravity  environment). 

In  addition  to  the  existing  FDA  approved  osteoporosis  drugs,  research  is 
ongoing  with  Osteoprotegerin,  another  naturally  occurring  protein  that  is  involved  in  the 
remodeling  process. 

c.  Other  Countermeasures 

Other,  more  exotic,  countermeasures  are  under  development  to  reduce  the 
effects  of  long-term  space  flight.  One  such  countermeasure  makes  use  of  vibrating 
plates.  Connecting  an  astronaut  to  a  plate  which  is  gently  vibrating  at  90  hz  for  10  to  20 
minutes  may  stimulate  osteoblasts  into  generating  new  bone.  Such  vibrations  have 
produced  near  normal  bone  formation  rates  in  studies  where  animals  were  prevented  from 
bearing  weight  in  certain  limbs.  The  control  animals  exhibited  a  92%  reduction  in  bone 
formation  when  not  exposed  to  the  vibrations  and  were  prevented  from  bearing  weight  in 
the  same  limbs.  Though  promising,  vibrating  plates  do  not  address  the  issues  of  muscle 
atrophy  and  degradation  of  aerobic  capacity. 20’2i 

Kathleen  M.  Ariss,  “Raloxifene  for  osteoporosis,”  WebMD.  2002,  Healthwise  Inc.,  14  Sept.  2002, 
http://mv.webmd.com/encvclopedia/article/1829.50935 

Ariss,  “Calcitonin  for  osteoporosis,”  WebMD.  2002,  Healthwise  Inc.,  14  Sept.  2002, 
http://mv.webmd.eom/encvclopedia/article/1829.50940 

!®  Ariss,  “Raloxifene  for  osteoporosis” 

Deglin,  Vallerand,  p.  23 

Ariss,  “Alendronate  for  osteoporosis,”  WebMD.  2002,  Healthwise  Inc.,  14  Sept.  2002, 
http  ://mv.  webmd.com/encvclopedia/article/ 1 829. 5093  7 

Ariss,  “Risendronate  for  osteoporosis,”  WebMD,  2002,  Healthwise  Inc.,  14  Sept.  2002, 
http://mv.webmd.com/encvclopedia/article/1829.50935 

20  Barry 

2!  “Astronaut  osteoporosis,”  BBCi,  15  May  2002,  British  Broadcasting  Corporation,  23  Sept.  2002, 
http://www.bbc.co.uk/science/tw/2002/mavl5osteoporosis.shtml 
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Another  countermeasure  that  is  under  development  is  known  as  Lower 
Body  Negative  Pressure  (LBNP).  An  LBNP  device  consists  of  a  partial  vacuum  chamber 
that  encompasses  the  legs,  feet,  and  pelvis,  which  seals  at  the  waist.  A  vacuum  of  30  to 
50  mmHg  is  pulled  in  the  chamber  resulting  in  a  footward  force,  which  through  design, 
can  be  the  equivalent  of  1  G.  Exercise  equipment,  such  as  a  treadmill,  can  be 
incorporated  within  and  around  the  LBNP  device  to  allow  an  astronaut  to  exercise  in  the 
equivalent  of  a  1  G  environment.  Exposing  the  body  to  a  vacuum  carries  with  it  the 
associated  risks  of  petechiae  (minute  hemorrhages  resulting  from  burst  capillaries22), 
hernia,  and  syncope  (brief  unconsciousness).  Eurthermore,  what  are  the  medical  risks  of 
long  term,  daily  exposure  of  the  lower  extremities  to  partial  vacuum? 

d.  Artificial  Gravity 

Using  a  centrifuge  to  provide  astronauts  with  artificial  gravity  in  space 
goes  to  the  root  of  the  problem,  weightlessness. 23  The  Centrifugal  Eorce  provided  by  a 
rotating  environment  simulates  the  presence  of  gravity.  This  rotating  environment 
provides  the  best  means  of  stimulating  the  natural  processes  that  maintain  the  body’s 
fluid  distribution,  muscle  tone,  and  bone  mass.  A  ‘normal’  exercise  routine  could  be 
enjoyed  on  ‘normal’  exercise  equipment.  A  normal  lifestyle  could  be  enjoyed  versus 
specialized  astronaut  food  and  sponge  baths. 

Additionally,  research  could  be  conducted  on  the  long-term  effects  of 
‘fractional-G’  on  the  human  body.  How  much  gravity  is  necessary  for  a  person  to 
maintain  physical  fitness  and  the  reformation  process?  Would  gravity  exposure  at 
scheduled  times  as  part  of  the  astronaut’s  daily  routine  be  sufficient  to  maintain  their 
health?  Outside  of  drop-tubes  and  parabolic  flight  paths  that  simulate  the  microgravity 
environment  for  a  few  seconds,  research  into  artificial  gravity  at  levels  of  centrifugal 
force  less  than  1  G  is  not  possible  on  Earth.  The  installation  of  a  centrifuge  aboard  a 
spacecraft  is  the  only  means  by  which  such  research  can  be  conducted. 


22  “Petechia,”  The  American  Heritage  Dictionary.  2““^  College  Ed.,  Boston:  Houghton-Mifflin  Co., 
1985,  p.  927 

23  Hullander,  Barry 
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B.  WHY  MODEL  A  SPACE  FLIGHT  CENTRIFUGE? 

Human  beings  are  not  designed  to  live  in  a  rotating  environment.  A  positive  side 
effect  of  a  rotating  environment  is  known  as  a  somatogravic  illusion:  a  false-perception 
of  attitude  due  to  prolonged  angular  motion.24  Through  this  illusion,  an  astronaut  will 
not  perceive  the  rotation  of  the  centrifuge,  only  the  force  of  its  artificial  gravity. 
However,  when  there  is  a  change  in  angular  velocity,  or  the  astronaut  changes  his  radial 
position,  a  “Coriolis  Illusion”  is  generated  that  contains  an  element  of  rotation  in  the 
plane  of  the  head  movement.^s  The  causes  and  effects  of  the  “Coriolis  Illusion”  will  be 
discussed  in  a  later  section. 

The  use  of  a  centrifuge  in  space  bears  associated  costs  and  risks.  More  to  the 
point,  can  a  centrifuge  be  constructed  and  put  in  space  at  a  cost  that  would  make  its 
realization  feasible?  The  engineering  behind  such  a  centrifuge  is  nothing  new,  the 
International  Space  Station  (ISS)  when  finished  will  have  a  centrifuge  as  part  of  the 
Centrifuge  Accommodation  Module  (CAM).  The  human  factor,  however,  is  not  so  well 
understood.  In  other  words,  what  is  the  best  balance  between  centrifuge  size  (cost)  and 
an  astronaut’s  ability  to  tolerate  the  rotational  environment  such  a  centrifuge  would 
impose?  Such  costs  and  benefits  must  be  weighed  not  only  in  dollars,  but  also  in  terms  of 
health,  time,  and  welfare  of  the  astronauts  it  would  benefit.  For  example,  suppose  that 
the  ISS  3  man  crew  were  able  to  gain  4.5  man-hours  every  day  by  reducing  aerobic 
exercise  requirements  from  2  hours  daily  to  30  minutes  daily^^  through  the  use  of  a 
centrifuge  to  maintain  their  health.  What  would  that  extra  time  be  worth? 

Construction  of  a  centrifuge  in  space  large  enough  for  human  occupation  is 
currently  in  the  realm  of  science  fiction.  The  cost  and  risk  inherent  to  space  flight  make 
such  an  endeavor  too  costly  to  consider  in  the  near-term.  However,  a  smaller  centrifuge, 
being  developed  for  the  purpose  of  performing  experiments  with  animal  subjects  is  being 


24  John  Ernsting,  A.N.  Nicholson,  D.J.  Rainford,  Aviation  Medicine.  3’‘*  Ed.,  Oxford:  Butterworth- 
Heinemann,  1999,  p.  438 

25  Ernsting,  Nicholson,  Rainford,  p.  434 

26  “ACSM  Guidelines  For  Healthy  Aerobic  Activity,”  1 1  Sept.  2000,  American  College  of  Sports 
Medicine.  American  College  of  Sports  Medicine,  14  Sept.  2002,  http://www.ascm.org/pdf/Guidelines.pdf 
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developed  for  the  ISS.  The  CAM,  whieh  will  be  small  enough  to  be  earned  aboard  the 
Spaee  Shuttle,  will  be  large  enough  to  allow  experiments  on  small  animals. 

The  usefulness  in  a  Centrifuge  Simulation  beeomes  apparent  when  one  eonsiders 
the  need  for  artifieial  gravity  during  long-term  spaee  flight,  the  magnitude  of  undertaking 
involved  in  building  a  human  sized  eentrifuge  in  spaee,  and  the  need  to  design 
experiments  for  the  ISS  Centrifuge  that  will  most  aeeurately  reproduee  the  desired  human 
experienee  for  its  test  subjeets. 

The  purpose  of  this  thesis  is  to  allow  the  user  to  eonduet  a  unified  simulation  of 
the  engineering  eost  of  a  eentrifuge  versus  an  astronaut’s  ability  to  tolerate  the  rotational 
environment  that  sueh  a  eentrifuge  would  impose. 
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II.  THESIS 


A.  HUMAN  DISORIENTATION  AS  A  FACTOR  IN  SPACECRAFT 

CENTRIFUGE  DESIGN 

1.  Why  Is  Human  Disorientation  Such  an  Important  Factor? 

Disorientation  is  manifested  when  the  Visual  and  Vestibular  System  reeeives  a 
provoeative  stimulus. 27  In  other  words,  the  eyes  are  reporting  something  different  than 
the  sense  of  orientation.  This  sensory  oonfliet  leads  to  feelings  that  range  from  general 
malaise,  disorientation,  nausea,  and  ean  eause  emesis  (vomiting28)  and  ineapaeitation. 

The  Vestibular  System  eonsists  of  the  Semi-Cireular  Canals  that  detect  changes  in 
angular  velocity,  and  the  Otoliths  that  detect  linear  acceleration.  These  provocative 
stimuli  can  be  specified  according  to  the  sensory  systems  that  are  involved.  A  visual- 
vestibular  conflict  occurs  when  the  eyes  and  the  vestibular  receptors  report  incompatible 
information.  An  intravestibular  conflict  occurs  when  there  is  a  mismatch  in  the 
information  reported  from  the  semi-circular  canals  and  the  Otoliths. 29 

The  physical  characteristics  of  the  semi-circular  canals  are  well  known  and  have 
been  derived  through  dynamic  analysis. The  three  canals  are  firmly  coupled  to  the 
skull  in  order  to  experience  the  same  accelerations  as  the  head.^i  As  shown  in  Figure  4, 
each  canal  contains  a  fluid,  endolymph,  which  rotates  within  the  canal  whenever  the  skull 
rotates  in  space.  The  endolymph,  flowing  through  a  smooth  bore,  enjoys  laminar  flow 
making  the  flow  resistance  linearly  dependent  on  velocity.  Movement  of  the  endolymph 
deflects  the  Cupula,  a  hair-cell  transducer,  within  the  ampulla.  The  Cupula  also  acts  as  a 
weak  spring  and  will  restore  itself  to  zero  deflection  in  the  angular  velocity  is  at  a  steady 
state.  Deflection  of  the  Cupula  produces  a  neural  signal  that  informs  the  brain  of  change 


27  Emsting,  Nicholson,  Rainford,  p.  459-461 

28  “Emesis,”  The  American  Heritage  Dictionary,  p.  448 

29  Emsting,  Nicholson,  Rainford,  p.  459-461 

20  Milsum,  p.  186 

21  Emsting,  Nicholson,  Rainford,  p.  434 
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in  the  person’s  spatial  orientation  32’33  mathematieal  properties  of  this  system  are 
laid  out  in  Appendix  C. 


'f'  (Head) 


Figure  4:  Diagram  Of  A  Semi-Circular  Canap4 


Of  partieular  interest  to  this  thesis  is  the  speeial  ease  of  a  Cross-Coupled,  or 
Coriolis,  Stimulation  where  the  individual  is  being  rotated  about  a  partieular  axis  and 
then  moves  the  head  in  a  manner  that  produees  a  ehange  in  angular  veloeity.  The  result 
of  Coriolis  Stimulation  is  to  produee  an  erroneous  signal  of  a  turn  about  an  axis  that  does 
not  aeeord  with  either  the  axis  of  rotation,  or  the  axis  of  movement.  This  erroneous 
signal  persists  after  the  movement  has  been  eompleted  due  to  the  faet  the  Cupula  requires 
10  seeonds  or  more  to  return  to  equilibrium.  During  this  time  the  Otoliths  sense  the 
eorreet  attitude  of  the  individual  with  relation  to  apparent  gravity.  This  mismateh 
produees  a  potent  stimulus  for  indueing  motion  siekness  to  whieh  all  individuals  with  an 
intaet  Vestibular  System  may  sueeumb  if  the  angular  veloeity  of  the  rotating  environment 
and  amplitude  of  head  movement  are  high  enough,  and  there  are  a  suffieient  number  of 

22  Milsum,  p.  186 

23  Emsting,  Nicholson,  Rainford,  p.  426 


24  John  H.  Milsum,  Biological  Control  System  Analysis.  New  York:  McGraw-Hill,  1966,  p.  186 
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repetitions. 35  Seasickness  is  a  classic  example  of  this  effect,  particularly  on  small  boats 
in  rough  seas. 

The  disorientation  resulting  from  radial  translation  in  a  rotating  environment 
could  potentially  cause  a  great  deal  of  disorientation  and  discomfort,  which  could  inhibit 
an  astronaut’s  ability  to  perform  normal  functions.  Movements  that  on  the  ground  that 
are  taken  for  granted  such  as  moving  from  a  sitting  to  a  standing  position,  walking, 
exercise,  and  climbing  a  ladder  would  become  extremely  burdensome  and  uncomfortable. 

2.  Modeling  Engineering  Cost  Versus  Coriolis 

Due  to  the  expense  of  placing  objects  in  orbit,  there  will  be  a  strong  temptation  to 
keep  the  size  of  a  centrifuge  for  use  by  astronauts  to  a  minimum.  However,  as  discussed 
in  the  previous  section,  Coriolis  Stimulus  can  have  a  profound  effect  on  an  individual. 
Therefore,  if  a  centrifuge  were  to  be  employed  in  a  long  term  mission  in  space,  the  effect 
of  Coriolis  must  be  balanced  with  cost  in  the  design  process  so  that  astronaut  will  be  able 
to  function  with  an  acceptable  level  of  discomfort. 

To  appropriately  model  Engineering  Cost  and  Coriolis  two  competing  models 
must  be  designed.  For  Engineering  Cost,  a  model  of  the  external  characteristics  of  a 
centrifuge,  namely  its  size  and  power  requirements,  must  be  designed.  The  mathematical 
rigor  behind  the  development  of  the  Centrifuge  Model  is  shown  in  Appendix  A.  For 
Coriolis,  a  model  of  the  internal  characteristics  of  a  centrifuge  and  the  semi-circular 
canals,  namely  Centrifugal  Force,  Coriolis  Force,  and  the  dynamic  response  of  the 
Cupula,  must  be  designed.  The  mathematical  rigor  behind  the  development  of  the 
Human  Disorientation  model  is  shown  in  Appendices  B  and  C. 
a.  Centrifuge  Model 

The  construction,  assembly,  and  flight  of  a  centrifuge  for  use  in  space  is 
not  a  trivial  exercise.  The  development  of  a  computer  model  that  considers  centrifuges  of 
various  sizes  must  take  into  account  the  effect  such  a  centrifuge  would  have  on  the 
spacecraft  of  which  it  would  be  a  vital  component.  However,  a  sizable  centrifuge  has  not 


35  Emsting,  Nicholson,  Rainford,  p.  461-462 
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yet  been  flown  in  space.  The  only  example  is  the  proposed  design  for  the  ISS  Centrifuge 
the  characteristics  of  which  are  shown  in  Table  1. 


Description 

Specification 

Location 

Centrifuge  Accomodation  Module 

Mass 

2700  kg 

Rotating  Radius 

1.25  m 

Habitat  Size 

19  in  Wx  28  in  Hx  24  in  D 

Habitat  Mass 

116  kg 

Habitat  Quantity 

2 

Apparent  Gravity 

0.01-2  G  (0.01  G  increments) 

Spin-Up  Time 

5-60  min 

Table  1:  International  Space  Station  Centrifuge  Specifications^^ 


“Centrifuge,”  2  July  2002,  Space  Station  Biological  Research  Project.  NASA  Ames  Research 
Center,  6  July  2002,  http://brp.arc.nasa.gov/GBL/centriTC.html 
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Although  the  ISS  Centrifuge  has  not  been  flown  yet,  it  does  provide  a 
baseline.  This  baseline  can  be  used  to  determine  the  effect  similar  centrifuges  of  various 
sizes  would  have  on  a  potential  spacecraft.  The  model  constructed  from  this  baseline 
would  assume  that  the  centrifuge  would  be  a  rotating  cylinder,  as  shown  in  Figure  1, 
designed  to  reproduce  the  dimensional,  mass,  apparent  gravity,  and  spin-up 
characteristics  of  the  ISS  Centrifuge.  This  model  would  be  scalable  to  test  the 
characteristics  of  cylindrical  centrifuges  of  various  sizes.  The  outputs  of  this  model  will 
be  mean  power  cost  and  graphs  of  angular  velocity,  current,  and  power  vs.  time. 

This  thesis  only  considers  cylindrical  centrifuges  rotating  as  shown  in 
Figure  5.  Some  futurists  and  visionaries  have  also  discussed  alternative  centrifuge 
architectures  such  as  radial-arm  and  tethered  centrifuges.  A  radial  arm  centrifuge  would 
consist  of  two  or  more  modules  mounted  at  the  end  of  rigid  spokes  connected  to  a  central 
hub  about  which  the  entire  apparatus  would  spin.  A  tethered  centrifuge  would  consist  of 
two  or  more  modules  connected  by  a  cable  that  would  be  held  rigid  by  the  centrifugal 
force  of  the  spiiming  bodies. These  alternate  centrifuge  architectures  present  different 
technical  characteristics  and  were  considered  to  be  beyond  the  scope  of  this  thesis. 

b.  Human  Disorientation  Model 

Consider  an  astronaut  within  a  centrifuge  operating  at  a  constant  angular 
velocity.  If  the  astronaut  is  stationary  relative  to  his  surroundings  he  feels  gravity  as  if  he 
was  on  Earth.  This  artificial  gravity  is  the  result  of  a  somatogravic  illusion  of  normal 
gravity  generated  within  a  centrifuge  because  the  apparent  gravitoinertial  force  vector  is 
normal  to  the  floor  on  which  he  is  standing,  as  in  Figure  6.  A  somatogravic  illusion  is  a 
phenomenon  of  the  Vestibular  System  that  produces  a  false  perception  of  actual 
gravitational  force  due  to  prolonged  exposure  to  an  atypical  force  vector.38  This  means 
that  the  astronaut  only  feels  the  centrifugal  force,  not  the  angular  velocity. 

As  shown  in  Equation  (B.14)  from  Appendix  B,  movement  by  the 
astronaut  in  the  axial  plane  (a  change  in  position  along  the  axis  of  rotation)  has  no  effect 
on  the  astronaut.  Movement  in  the  angular  plane  (a  change  in  position  around  the  axis  of 

Theodore  W.  Hall,  “After  Skylab,  1973-1991,”  The  Arehiteeture  of  Artifieial-Gravity 
Environmentsfor  Long-Duration  Spaee  Habitation.  Ann  Arbor:  University  of  Miehigan,  1994,  24  Sept. 
2002,  http://wwwO.areh.cuhk. edu.hk/~hall/ag/Dissertation/l  4.htni 

3^  Ernsting,  Nicholson,  Rainford,  p.  437 
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rotation)  produces  a  Coriolis  Force  that  has  the  effect  of  increasing  or  decreasing  the 
apparent  gravitational  force.  Movement  in  the  radial  plane  (change  in  distance  from  the 
axis  of  rotation)  produces  a  Coriolis  Force  that  has  the  effect  of  both  a  shear  force  in  the 
angular  plane  and  a  torque  in  the  axial  plane.  Previous  studies  conducted  where  human 
test  subjects  were  placed  into  a  slowly  rotating  room  underscored  the  detrimental  effects 
of  translation  in  the  radial  plane  of  a  rotating  environment.39 


Figure  6:  An  Astronaut’s  Perception  Of  The  Rotating  Environment 

When  the  astronaut  moves  in  the  radial  plane,  he  receives  a  Coriolis 
Stimulation  and  detects  a  change  in  angular  velocity  with  his  semi-circular  canals.'^*’  The 
two  sets  of  semi-circular  canals  reside  on  each  side  of  the  head  in  the  inner  ear  (illustrated 
in  Figure  7  and  Figure  8).  Each  semi-circular  canal  is  capable  of  detecting  angular 
acceleration  in  three  dimensions  via  the  anterior  vertical,  posterior  vertical,  and  lateral 
canals;  labeled  “pv,”  “av,”  and  “1”  in  Figure  8.  The  brain  processes  the  information  from 
each  semi-circular  canal  and  is  able  to  accurately  sense  the  plane,  direction,  and 
magnitude  of  any  change  in  angular  motion  within  certain  limits.  Although  the  semi¬ 
circular  canals  are  not  aligned  with  the  pitch,  roll,  and  yaw  axes  the  brain  processes  the 
information  from  both  canals  as  a  unified  model  of  the  body’s  angular  motion. 

39  Percival  McCormack,  Personal  Interview,  7  Feb.  2002 
Ernsting,  Nicholson,  Rainford,  p.  425 
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Therefore,  the  alignment  of  the  semi-cireular  canals  is  trivial  and  only  head  alignment 
should  be  considered.'^! 


Figure  7:  The  Inner  Ear'^^ 


Figure  8:  Semi-Circular  Canal 
Orientation'*^ 


As  shown  in  equation  (B.14),  in  Appendix  B,  Coriolis  Stimulation  due  to 
radial  translation  only  occurs  in  one  plane.  Therefore,  only  one  semi-circular  canal  needs 
to  be  modeled.  The  effects  of  linear  acceleration  due  to  Coriolis  are  not  considered  in 
this  paper  because  they  are  assumed  to  be  predictable  and  astronauts  should  be  able  to 
adapt  to  them  easily. 


B,  COMPUTER  MODELS 

Placing  objects  in  orbit  is  extremely  expensive.  As  a  result  there  will  be  a  strong 
temptation  to  keep  the  size  of  any  centrifuge  considered  for  use  in  space  by  astronauts  to 
a  minimum  without  consideration  of  the  profound  effect  Coriolis  Stimulation  can  have  on 
an  individual.  In  an  effort  to  promote  the  use  of  a  centrifuge  in  long-term  missions  in 
space,  two  competing  models  were  designed,  appropriately  modeling  Engineering  Cost 
and  the  effect  of  Coriolis.  In  order  for  the  use  of  a  centrifuge  in  space  to  be  effective,  a 

4!  Emsting,  Nicholson,  Rainford,  p.  426 

42  Emsting,  Nicholson,  Rainford,  p.  424 

43  Emsting,  Nicholson,  Rainford,  p.  426 
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balance  must  be  struek  between  its  Engineering  Cost  and  the  effeet  of  Coriolis  on  its 
occupants  must  be  balaneed  in  the  design  proeess  so  that  an  astronaut  will  be  able  to 
function  with  an  acceptable  level  of  diseomfort. 

1.  Building  a  Computer  Model  Using  LabView 

National  Instruments  LabView  6.0  was  seleeted  as  the  program  within  whieh  the 
necessary  models  were  to  be  developed.  Its  advantages  were  its  graphieal  interface,  its 
stability,  the  modularity  of  its  eomponents,  ease  in  debugging,  the  portability  of  the 
programs  developed  using  it  to  multiple  platforms,  and  its  ability  to  run  on  the  several 
platforms  of  various  eapabilities. 

2.  Dual  Analysis 

The  equations  of  motion,  as  developed  in  the  appendiees,  were  solved  using  two 
different  mathematieal  approaehes.  First,  a  eheek  ease  was  developed  using  a  step  inputs 
and  Laplaee  Transforms  to  aehieve  an  analytieal  solution  of  the  model  as  a  function  of 
time.  Seeond,  eaeh  model  was  written  in  state-variable  form,  diseretized  and 
programmed  in  LabView  to  allow  numerieal  solutions  for  arbitrary  input  eonditions. 
Eaeh  analysis  was  eompared  using  separate  eomputer  programs  developed  in  LabView  to 
determine  if  eaeh  analysis  was  produeing  the  same  numerieal  results. 

3.  Experimental  Validation 

A  simple  validation  experiment  was  eonducted  to  get  a  ‘seat-of-the-pants’  feel  for 
the  results  of  the  Human  Disorientation  Model.  Using  a  Merry-Go-Round  was  the 
simplest  means  to  determine  how  radial  movement  of  the  head  feels  in  a  rotating 
environment  roughly  the  size  of  the  ISS  Centrifuge.  A  7  ft.  6.5  in.  Merry-Go-Round, 
pietured  in  Figure  9,  was  found  on  the  Former  Fort  Ord  in  Marina,  CA  near  the  Naval 
Postgraduate  Sehool;  the  outer  radius  of  whieh  was  within  10  em  of  the  rotating  radius  of 
the  ISS  Centrifuge.  The  experiment  was  designed  to  be  run  with  a  small  number  of 
people  of  various  military  baekgrounds  who  were  eolleagues  of  the  author  in  the  Spaee 
Systems  Operations  Currieulum  at  the  Naval  Postgraduate  Sehool.  The  experiment  was 
also  designed  to  be  run  at  levels  of  Centrifugal  Foree  that  eould  be  withstood  by  a  person 
sitting  on  a  Merry-Go-Round. 
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Experiment  Run 

Centrifugal  G 

10  RPM 

0.13 

20  RPM 

0.52 

Table  2:  Centrifugal  Force  Produced  During  Each  Run 

The  experiment  consisted  of  two  runs  per  test  subject.  The  first  run  would  be 
conducted  at  10  RPM,  and  the  second  at  20  RPM.  The  angular  velocities  were  selected 
to  keep  the  centrifugal  force  at  a  level  that  could  be  managed  by  an  individual  sitting  on 
an  open  platform  Merry-Go-Round  as  shown  in  Table  2. 


Figure  9:  The  Apparatus  Used  In  The  Validation  Experiment 


Figure  10:  Test  Subject’s  Perspective 


17 


Each  run  was  then  divided  into  three  phases  for  the  purposes  of  reeording  the 
amount  of  disorientation.  The  first  phase  eonsisted  of  the  spin-up  aeeeleration  and  first 
60  seeonds  of  steady  state  angular  veloeity.  The  seeond  phase  eonsisted  of  10  seeonds  of 
eyolieal  head  movement  at  about  0.5  Hz  through  10  em.  The  third  phase  eonsisted  of  the 
final  60  seeonds  and  the  spin-down  aeeeleration.  Additionally,  the  test  subjeets  were 
instrueted  to  stare  at  the  yellow  surveyor’s  flag  direetly  aeross  from  them  so  that  any 
oeulogyral  illusions  around  that  fixed  point  eould  be  reeorded  as  shown  in  Figure  10.  An 
arbitrary  seale  of  disorientation,  shown  in  Table  3,  was  devised  to  quantify  the 
disorientation  felt  by  the  test  subjeets.  Time  was  kept  using  a  stopwateh  and  angular 
veloeity  was  maintained  by  eounting  the  seeonds  on  the  stopwateh  between  passes  of  the 
yellow  surveyor’s  flag  shown  in  Figure  9. 


Disorientation  Fevel 

Assoeiated  Feeling 

1 

Normal 

2 

Slightly  Dizzy 

3 

Dizzy 

4 

Nauseated 

5 

Imminent  Emesis 

Table  3:  Qualitative  Disorientation  Scale 


4.  Expected  Results 

The  eomputer  models,  eombined  with  the  results  of  the  validation  experiment 
should  show:  (1)  that  disorientation  is  felt  when  the  head  is  moved  in  the  radial  plane  of  a 
rotating  environment,  (2)  that  the  disorientation  felt  is  related  to  a  feeling  of  twisting  in 
the  axial  plane  of  a  rotating  environment,  and  (3)  that  the  feeling  of  disorientation  is 
direetly  related  to  eentrifuge  radius  and  angular  veloeity.  The  validation  experiment 
allows  for  a  ealibration  of  predieted  Cupula  Veloeity  to  be  related  to  an  arbitrary  level  of 
diseomfort.  Additionally,  this  thesis  should  show  that  a  eentrifuge  for  providing 
astronauts  artifieial  gravity  without  signifieant  diseomfort  due  to  disorientation  is  feasible 
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and  could  be  a  future  upgrade  to  the  International  Spaee  Station,  part  of  a  future  space 
station,  or  a  component  of  a  spaceeraft  for  a  manned  interplanetary  flight. 
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III.  COMPUTER  MODEL 


A,  COMBINED  SIMULATION 

The  ability  to  develop  a  dynamic  model  of  a  centrifuge  and  its  effect  on  human 
disorientation  resides  in  the  mathematical  understanding  of  the  two  systems.  Because 
this  thesis  seeks  to  produce  a  unified  analysis  of  two  different  systems;  a  user  interface 
was  developed  that  allows  the  simultaneous  input  of  parameters  and  display  of  results. 
This  overview  of  the  two  models  is  called  the  Centrifuge  Front  Panel. 


/fftOIITV 

\PflllELA 

Centrifuge  Front  Panel.vi 

Centrifuge  State  Space. vi 

Astmi 

kStfiCl 

Human  Disorientation  State  Space. vi 

Table  4:  Major  Components  Of  The  Centrifuge  Front  Panel 

The  Centrifuge  Front  Panel  is  a  computer  program  written  using  LabVIEW 
This  “virtual  instrument,”  or  VI  ,  as  shown  in  Table  4,  incorporates  the  Centrifuge 
Dynamic  Model,  the  Human  Disorientation  Model,  and  their  various  sub-components  to 
display  the  engineering  costs  of  an  arbitrary  centrifuge  together  with  the  physical 
response  of  a  human  semi-circular  canal.  The  block  diagram  shown  in  Figure  1 1  details 
the  relationships  between  the  various  Vis  utilized  in  this  thesis. 

1,  Combined  Simulation  User  Interface 

The  user  interface  (known  as  the  front  panel)  of  the  combined  simulation,  shown 
below  in  Figure  12,  is  designed  to  allow  the  user  to  input  the  specific  parameters  of  a 
centrifuge  and  the  person  occupying  it  via  controls  with  the  results  displayed  on 
indicators.  The  individual  parameters  (known  as  controls)  are  available  for  user 
manipulation  when  up  and  down  arrows  are  depicted  next  to  the  numerical  readout. 
These  parameters  are  grouped  into  ‘clusters’  of  related  data  such  as  Simulation 
Parameters,  Electric  Parameters,  Centrifuge  Dimensions,  and  Human  Parameters.  The 

Lab  View.  CD-ROM,  Ver.  6i,  Austin  Texas:  National  Instruments  Corp.,  2000 
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output  paramters,  known  as  indicators,  of  the  simulation  took  the  form  of  plots  or 
numerieal  readouts.  Every  control,  cluster,  or  indicator  is  represented  on  a  diagram, 
which  is  where  LabVIEW programs  are  written. 


Frequency  & 
Omega  Converter 


Disorlenlatior 
State  Space 

T 

**\ 

L 

Semicircular 
Canal  Dyrramics 

Har>d| 

Centrifuge  Steady 

LabVtew 

Total  Moment  of 

State 

mean.vi 

Inertia 

yjjjU 
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Converter 

Disc  Mass  & 

HolkAv  Cylinder 

Rod  Mass  & 

Moment  of  Inertia 

Mass  &  MOI 

Moment  of  Inertia 

Figure  11:  Block  Diagram  of  the  Centrifuge  Front  Panel 


This  VI  is  a  combined  simulation  of 
Centrifuge  Dynamics  and  Human  Disorientation 
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Figure  12:  Centrifuge  Front  Panel 
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2,  Diagram 

a.  Centrifuge  Dynamic  Model 

The  diagram  of  the  Centrifuge  Front  Panel  eontains  a  Sequenee  Strueture 
that  consists  of  two  frames.  The  first  frame,  pictured  below  in  Figure  13,  collects  the 
user’s  inputs  for  the  Centrifuge  Dynamics  Model.  Those  inputs  are  the  Electric 
Parameters,  Centrifuge  Dimensions,  and  Simulation  Parameters  clusters.  Wiring  the 
clusters  from  the  front  panel  of  this  VI  into  a  linked  VI,  in  this  case  Centrifuge  State 
Space,  has  the  effect  of  reproducing  the  user’s  inputs  in  similarly  structured  control 
clusters  in  the  linked  VI.  One  output.  Angular  Velocity  (a  one-dimensional  array  indexed 
with  respect  to  time)  is  wired  to  a  Sequence  Local  that  passes  that  data  to  any  subsequent 
frames  for  later  use.  The  remaining  outputs;  scalars  representing  Centrifuge  Gravity  and 
Mass  and  vectors  representing  Power  vs.  Time  and  Angular  Velocity  vs.  Time  are 
connected  to  their  corresponding  display  indicators. 


Figure  13:  Front  Panel  Interaction  With  The  Centrifuge  Dynamic  Model 


b.  Human  Disorientation  Model 

The  second  frame,  pictured  in  Figure  14,  collects  the  inputs  for  the  Human 


Disorientation  Model.  Those  inputs  are  the  Human-Centrifuge  Interaction,  SCC 
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Characteristics,  Simulation  Parameters,  and  Human  Parameters  clusters  together  with  the 
Radial  Translation  scalar  and  the  Angular  Velocity  array  from  the  Centrifuge  Dynamics 
Model.  The  inputs  are  wired  to  their  counterpart  controls  in  the  linked  VI,  in  this  case 
Human  Disorientation  State  Space.  The  outputs  are  vectors  representing  Radial  Position 
vs.  Time,  Cupula  Velocity  vs.  Time,  and  Gravity  At  The  Head  vs.  Time  are  connected  to 
their  corresponding  display  indicators. 


y^fr^ 

Centrifuge  State  Space. vi 

If 

us^ 

Centrifuge  Steady  State.vi 

Freq  and  Omega  2  way  Converter. vi 

Mx 

Mean 

Mean.vi 

©•j 

Total  Moment  of  Inertia,  vi 

Table  5:  Major  Components  Of  The  Centrifuge  Dynamic  Model 
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The  details  of  the  construction  of  the  Centrifuge  Dynamic  Model  and  the 
Human  Disorientation  Model  Vis  will  be  detailed  in  the  proceeding  sections. 


B.  CENTRIFUGE  DYNAMICS  MODEL 

Because  the  Vis  listed  in  Table  5  deal  specifically  with  the  dynamics  of  a 
centrifuge,  more  information  can  be  displayed  in  the  user  interface.  The  Control  Clusters 
of  this  VI  are  structured  so  that  the  model  can  be  run  independently,  or  as  component  of 
the  combined  simulation. 

1,  Front  Panel 


This  VI  computes  the  dynamic  response 
of  a  centrifuge  based  on  an  arbitrary  input 
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Figure  15:  Centrifuge  Dynamic  Model  User  Interface 


The  front  panel  of  the  Centrifuge  Dynamics  Model,  shown  in  Figure  15,  is 
structured  so  that  the  user  can  not  only  control  the  simulation  (as  when  running  the 
combined  simulation),  but  can  also  vary  the  start-up  dynamics  so  that  different  means  of 

reducing  the  power  required  to  start  the  centrifuge  can  be  examined.  The  different  start- 
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up  functions  are:  Unit-Step,  Ramp,  Two-Step,  and  Three-Step.  Additional  indieators  are 
utilized  to  show  speeific  faeets  of  the  eentrifuge’s  engineering  cost,  specifically:  Mean 
Power  Cost,  Voltage  vs.  Time,  and  Current  vs.  Time. 

2,  Diagram 

a.  Centrifuge  Natural  Response  Characteristics 

The  diagram  of  the  Centrifuge  Dynamie  Model  eontains  a  Sequenee 
Strueture  that  eonsists  of  four  frames.  The  first  frame,  pietured  in  Figure  16,  eollects  the 
user’s  inputs  and  determines  the  Centrifuge  Steady  State  Response.  Those  inputs  are  the 
Eleetrie  Parameters,  Centrifuge  Dimensions,  and  Simulation  Parameters  elusters  ean  be 
entered  direetly  by  a  user  or  linked  from  a  superior  VI.  The  Centrifuge  Dimension 
cluster  of  this  VI  is  wired  into  the  linked  VI,  Total  Moment  of  Inertia  and  the  Formula 
Node  for  determining  Steady  State  Centrifugal  Foree.  The  Eleetrie  Parameters  cluster, 
together  with  the  output  of  the  Total  Moment  of  Inertia  VI,  is  rebundled  into  the  Physieal 
Parameters  eontrol  eluster  of  the  linked  Centrifuge  Steady  State  Response  VI.  The  output 
of  the  Centrifuge  Steady  State  Response  VI  is  rebundled,  together  with  the  Total  Mass 
Moment  of  Inertia,  into  the  Response  Characteristies  indieator  eluster.  Steady  State 
Angular  Veloeity  is  separately  unbundled  from  the  Response  Characteristies  eluster  to 
provide  an  input  to  the  Formula  Node  for  eomputing  Steady  State  Centrifugal  Foree. 
Total  Moment  of  Inertia  is  wired  to  a  Sequenee  Foeal,  and  the  Formula  Node  returns 
Steady  State  Centrifugal  Foree  for  its  eorresponding  indieator  on  the  Front  Panel. 

The  following  formula,  based  on  Equation  (B.18)  in  Appendix  B,  is  used 
to  eompute  the  Steady  State  Centrifugal  Force: 


go 
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b.  Unit  Step  Function  Start-Up  Response 


Logic 

Case  Response 

Input  Stop  Time  >  Current  Time 

Input  Start  Time  >  Current  Time 

False 

Input  Stop  Time  >  Current  Time 

Input  Start  Time  <=  Current  Time 

True 

Input  Stop  Time  <=  Current  Time 

Input  Start  Time  <=  Current  Time 

False 

Table  6:  Start-Up  Voltage  Case  Logic 


The  second  frame,  pictured  in  Figure  17,  collects  the  user’s  inputs  and 
computes  voltage  with  respect  to  the  start-up  function.  Those  inputs  are  the  Electric 
Parameters  and  Simulation  Parameters  clusters.  The  start-up  function  is  determined 

through  the  use  of  a  Ring  Control.  The  first,  and  default,  selection  of  the  Ring  Control  is 
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the  Unit  Step  Function.  The  Ring  Control  is  not  designed  to  be  manipulated  from  a 
superior  VI,  therefore  when  the  Centrifuge  Dynamic  Model  is  linked  within  another  VI,  it 
will  generate  its  responses  using  the  Unit  Step  Function  which  means  that  the  voltage 
instantaneously  changes  from  zero  to  its  steady  state  value. 


Figure  17:  Centrifuge  Dynamic  Start-Up  Response  to  a  Unit-Step  Function 


For  each  case,  a  simple  logic  tree,  shown  in  Table  6,  and  Case  Structure  is 
used  to  determine  when  to  apply  the  designated  start  function.  Voltage,  Time,  and  a 
Voltage  Applied  Boolean  are  indexed  into  one-dimensional  arrays  and  connected  to 
Sequence  Locals  for  later  use. 
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Figure  18:  Centrifuge  Dynamic  Start-Up  Response  To  A  Ramp  Function 


c.  Ramp  Function  Start-Up  Response 

The  second  case,  pictured  in  Figure  18,  uses  a  Ramp  Function  to  compute 


the  start-up  voltage.  The  Ramp  Function  utilized  is  shown  below: 


V  =V 

Start— Up  final 


^function  J 


d.  Two-Step  Function  Start-Up  Response 

The  third  case,  pictured  in  Figure  19,  uses  a  Two-Step  Function  to 
compute  the  start-up  voltage.  This  case  assumes  some  electronic  means  is  used  to  halve 
the  voltage  for  the  duration  of  the  start-up  function.  The  Two-Step  Function  utilized  is 
shown  below: 

F  =iF 

Start-up  r\  fino.1 
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e.  Three-Step  Function  Start-Up  Response 

The  fourth  case,  pictured  in  Figure  20,  uses  a  Three-Step  Function  to 
compute  the  start-up  voltage.  This  case  assumes  some  electronic  means  is  used  to  step 
the  voltage  by  thirds  for  the  duration  of  the  start-up  function.  The  first  half  of  the  Three- 
Step  Function  utilized  is  shown  below: 

F  =1f 

Start-up  ^  final 


The  second  half  of  the  Three-Step  Function  is: 

V  =-V 

start-up  ^  final 
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Figure  20:  Centrifuge  Dynamic  Start-Up  Response  to  a  Three-Step  Function 


f  Initial  Prediction  of  Current  and  Angular  Velocity 

The  third  frame,  pietured  in  Figure  21,  eonsists  of  a  nested  Sequenee 
Structure  that  has  three  frames.  The  first  nested  frame  takes  user  inputs  and  the  current 
state  of  angular  velocity  and  current  to  compute  the  change  in  current  and  angular 
velocity  and  predict  the  new  state  in  angular  velocity  and  current.  The  user  inputs  are  the 
Electric  Parameters  and  Simulation  Parameters  clusters.  Moment  of  Inertia,  Voltage, 
Time,  and  the  Voltage  Applied  Boolean  are  brought  forward  into  this  frame  using 
Sequence  Locals.  The  Voltage  Applied  Boolean  is  used  to  determine  whether  or  not 
voltage  is  being  applied. 

The  formulae  used  in  Figure  21  are  based  on  the  State  Space  Analysis 
found  in  Equation  (A.  123)  in  Appendix  A.  Change  in  current  is  predicted  by: 

?  R  .  V 

i  =  — i  — —0)-\ — 

ILL 
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Change  in  angular  velocity  is  predicted  if  voltage  is  applied  by: 

^  k,  .  B 

0)  =  —l — 0) 

J  J 

Otherwise,  the  change  in  angular  velocity  without  voltage  is: 

-  B 
(0  =  — (O 

J 


Figure  21:  State  Space  Prediction  of  Current  and  Angular  Velocity 


The  predicted  current  when  voltage  is  applied: 

i  =/  +  /At 

Otherwise,  without  voltage  the  predicted  current  is  zero.  Predicted 
angular  velocity  is  determined  by: 

a)  =  a)+d)dd 

'  "  .V 

The  results  of  the  above  formulae:  z  ,  z ,  z»,  and  Q)  are  passed  to  the  next 
frame  using  Sequence  Locals. 
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g.  Correction  of  Predicted  Current  and  Angular  Velocity 
The  second  nested  frame  takes  user  inputs  and  the  predicted  state  of 
angular  velocity  and  current  to  correct  the  predicted  state  of  change  in  current  and 
angular  velocity  and  produce  the  new  state  of  angular  velocity  and  current.  Again,  the 
user  inputs  are  the  Electric  Parameters  and  Simulation  Parameters  clusters.  Moment  of 
Inertia,  Voltage,  Time,  and  the  Voltage  Applied  Boolean  are  brought  forward  into  this 
frame  using  Sequence  Locals.  Additionally,  the  predicted  states  are  passed  into  this 
frame  using  Sequence  Locals.  The  Voltage  Applied  Boolean  is  used  to  determine 
whether  or  not  voltage  is  being  applied. 


Figure  22:  Correction  of  Predicted  State  to  Produce  New  State 


The  formulae  used  in  Ligure  22  are  based  on  those  found  in  Equation 
(A.  123)  in  Appendix  A.  Change  in  current  is; 


R,  k.  V 

■■ — l - -CO+  — 

ILL 
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Change  in  angular  velocity  is  when  voltage  is  applied; 

.  k,,B. 

CO  =  —l - CO 

J  J 

Otherwise,  the  change  in  angular  velocity  without  voltage  is: 

B  , 

0)  =  — CO 

J 


The  corrected  prediction  of  current  when  voltage  is  applied: 

At 


=h  + 


f  ?  o 

I  +l 


V  2  , 


Otherwise,  without  voltage  the  predicted  current  is  zero.  The  corrected 
prediction  of  angular  velocity  is  determined  by: 


6)  +  6)^ 

V  ^  y 


At 


In  addition  to  the  current  and  angular  velocity  predictions,  instantaneous 
power  is  computed  in  this  frame: 


P  =  i  V 


The  results  of  the  computations  conducted  within  this  frame:  Power, 
Current,  and  Angular  Velocity  are  wired  outside  of  the  nested  Sequence  Structure, 
indexed  versus  time  into  one-dimensional  arrays,  and  connected  to  Sequence  Locals  for 
use  in  a  later  frame. 


h.  Display  of  Results 

In  the  final  frame  of  the  Centrifuge  Dynamic  Model,  the  results  of  the 
State  Space  Analysis:  Power,  Current,  and  Angular  Velocity  are  bundled  with  a  one¬ 
dimensional  Time  array  and  plotted  on  the  Front  Panel.  Additionally,  Mean  Power  Cost 
is  computed  using  mean.vp^  and  displayed  as  a  scalar  together  with  Outer  Radius.  A 
one-dimensional  of  the  Angular  Velocity  results  is  connected  to  a  hidden  indicator  for 
linking  to  the  Human  Disorientation  Model. 


45  “mean.vi”,  CD-ROM,  Lab  VIEW  6i.  Austin  Texas:  National  Instruments  Corp.,  2000 
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C.  CENTRIFUGE  STEADY  STATE  RESPONSE 
1,  Front  Panel 


Centrifuge  Steady  State. vi 

f^<«) 

Freq  and  Omega  2  way  Converter. vi 

Degrees  and  Radians  2  way  Converter.vi 

Table  7:  Major  Components  of  the  Centrifuge  Steady  State  Model 

The  front  panel  of  the  Centrifuge  Steady  State  Response,  shown  in  Figure  24,  is 
not  designed  for  direct  user  manipulation.  Rather,  this  VI  is  designed  to  be  a  subordinate 
VI  nested  within  the  Centrifuge  State  Space  VI.  The  variables  within  the  input  Control 
Cluster,  Physical  Parameters,  are  linked  from  the  Centrifuge  Dynamic  Model  providing 
the  information  necessary  to  compute  the  steady  state  characteristics  of  a  centrifuge  in 
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motion.  The  output  Indicator  Clusters:  Steady  State  Parameters  and  Response 
Characteristics  link  the  results  of  this  VI  to  the  Centrifuge  Dynamic  Model.  The  major 
components  of  this  VI  are  listed  in  Table  7. 


This  VI  Computes  the  steady  state  response  and 

damping  chacteristics  of  the  system 

Physical  Parameters 

Steady  State  Parameters 

Bus  Voltage,  V 

Steady  State  Angular  Velocity,  deg/sec 

28.000 

15.884 

Steady  State  Angular  Velocity,  rad/sec 

Friction  Losses,  N 

0.277 

1.000 

Steady  State  Angular  Velocity,  hz 

0.044 

Motor  Resistance,  ohms 

Steady  State  Angular  Velocity,  RPM 

y.  100.000 

2.647 

Steady  State  Current,  A 

Motor  Inductance,  H 

0.028 

10.000 

Response  Characteristics 

Mass  Moment  of  Inertia,  kg^m? 

Steady  State  Angular  Velocity,  rad/sec 

100.000 

0.277 

Steady  State  Current,  A 

T  orque  Proportionality  Constant,  k_t 

0.028 

rjq.ooo 

Damping  Ratio 

0.173 

EMF  Proportionality  Constant,  k_f 

rj  1.000 

Natural  Frequency,  rad/sec 

31.780 

Figure  24:  Centrifuge  Steady  State  Model  User  Interface 


2,  Diagram 

a.  Steady  State  Angular  Velocity 

The  diagram  of  the  Centrifuge  Steady  State  Response  contains  a  Sequence 
Structure  that  consists  of  five  frames.  The  first  frame,  pictured  in  Figure  25,  collects  the 
linked  inputs  and  determines  the  Centrifuge’s  Steady  State  Angular  Velocity.  The  Steady 
State  Angular  Velocity  computed  in  this  frame  is  wired  into  a  Sequence  Local  for  use  in 
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later  frames.  The  fifth  frame,  not  pictured,  converts  Steady  State  Angular  Velocity  into 
four  commonly  used  units  (rad/sec,  deg/sec,  hz,  and  RPM)  all  of  which  are  bundled  for 
display  in  the  Steady  State  Parameters  Indicator  Cluster.  The  formula  used  in  this  frame 
is  based  on  Equation  (A.  1 5)  in  Appendix  A; 

[RB  +  k^kf) 


iPhj^al  Parameters] 


— iFrictionI - 

]Resi$lance1 


This  formula  node  computes  steady  state  angular  velocity  in  rad/sec 
den  =  ( R  “  B  )  +  ( k_t  -  k_f ); 


onnega_ss  =  ( V  “  k_l  )/den. 


[  /?  5  +  A:,  k,  ] 

:  sL 


ISlead^)  State  Parameters! 


□ 

□ 

□ 

□ 

n 

Natural  Frequeni:y 
rad/sec 

Steady  State  Angular  Velocity 
rad/sec 

Response  Characteristics! 


current_$$  h 


-[Natural  Frequency,  rad/seci — 
- IDampinq  Ratio! - 


Figure  25:  Determination  Of  Centrifuge  Steady  State  Angular  Velocity 


b.  Steady  State  Current 

The  second  frame,  pictured  in  Figure  26,  collects  the  linked  inputs  and 
determines  the  Centrifuge’s  Steady  State  Current.  The  Steady  State  Current  computed  in 
this  frame  is  wired  into  a  bundle  for  display  in  the  Steady  State  Parameters  Indicator 
Cluster.  The  formula  used  in  this  frame  is  based  on  Equation  (A.  17)  in  Appendix  A; 
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L,  = 


B 


[RB  +  k,kf) 


V 


Figure  26:  Determination  of  Steady  State  Current 


c.  Natural  Frequency 

The  third  frame,  pictured  in  Figure  27,  collects  the  linked  inputs  and 
determines  the  Centrifuge’s  Natural  Frequency.  The  Natural  Frequency  computed  in  this 
frame  is  wired  into  a  bundle  for  display  in  the  Response  Characteristics  Indicator  Cluster. 
The  formula  used  in  this  frame  is  based  on  Equation  (A.25)  in  Appendix  A; 


(0„  = 


^  BR  +  k^kf'^ 


LJ 
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Figure  27:  Determination  of  Natural  Frequency 


d.  Damping  Ratio 

The  fourth  frame,  pictured  in  Figure  28,  collects  the  linked  inputs  and 
determines  the  Centrifuge’s  Damping  Ratio.  The  Damping  Ratio  computed  in  this  frame 
is  wired  into  a  bundle  for  display  in  the  Response  Characteristics  Indicator  Cluster.  The 
formula  used  in  this  frame  is  based  on  Equation  (A.25)  in  Appendix  A; 


(o„  = 


BR  +  k^kf'^ 
LJ 
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Figure  28:  Determination  of  Damping  Ratio 


D.  HUMAN  DISORIENTATION  MODEL 


/ftSTftTE 

Kstfici 

Human  Disorientation  State  Space. vi 

Freq  and  Omega  2  way  Converter. vi 

see 

Hand 

see  Hand  Math.vi 

Table  8:  Major  Components  of  The  Human  Disorientation  Model 

As  with  the  Centrifuge  Dynamic  model,  the  Vis  listed  in  Table  8  deal 
specifically  with  the  dynamics  of  a  Semi-Circular  Canal,  more  information  can  be 
displayed  in  the  user  interface.  The  Control  Clusters  of  this  VI  are  structured  so  that  the 
model  can  be  run  independently,  or  as  component  of  the  combined  simulation. 
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1. 


Front  Panel 


This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 


see  eharacteristics  Human-CGntrifugG  Interaction  eurrentTime  Point 


Restoring  Force,  k/J,  1/sec 
rJiTDOd 

Friction  Losses,  b/J,  1  /sec 
^200.00 

Initial  Radial  Position,  m 
;j(230 

Function  Speed,  Hz 
^0.50 

12000 

Number  Of  Points 

12001 

Current  Time,  sec 

Simulation  Parameters 

Human  Parameters 

120.00 

Input  Start  Time,  sec 

Function  Type 

Angular  Velocify,  rad/sec 

;jSland8;Sil  1 

Input  Start  Function  Duration,  sec  Subi.  Function  Start  Time,  sec 

1.0472 

rJA.QO 

Sim  Stop  Time,  sec 

;j|B0,00 

Subiect  Function  Duration,  sec 

Independent  — 

'.ymm 

i 

Sim  Sample  Interval,  sec 

^JIo.oiooF 

Input  Stop  Time,  sec 

;j|  20.00 

Radial  Translation,  m 

j  0.1000 

Angular  Velocity,  RPM 
:)  10.0000 

Integrated 

eoriolis,  N 
40.0-1 
20.0-! 
I  0-0- 
1  -20.0-1 


.JAiliwAili 

iniiiiniil 
iiiiiinnii 

, _ mmimimir 
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Time 
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75.000  80.000  85.000  90.000 


Figure  29:  Human  Disorientation  Model  User  Interface 


The  front  panel  of  the  Human  Disorientation  Model,  shown  in  Figure  29,  is 
structured  so  that  the  user  can  not  only  control  the  simulation  (as  when  running  the 
combined  simulation),  but  can  also  employ  various  radial  motion  functions  of  the  test 
subject  so  that  different  types  of  motion  can  be  tested.  The  motion  functions  are;  At- 
Rest,  Stand  &  Sit,  Up  Ladder,  and  Down  Ladder.  Additional  indicators  are  utilized  to 
show  specific  facets  of  the  human  disorientation,  specifically:  Cupula  Deflection  vs. 
Time,  Radial  Velocity  vs.  Time,  and  Coriolis  vs.  Time.  Finally,  a  toggle  switch  is 
employed  to  allow  the  user  to  change  the  model  between  independent  and  integrated 
modes.  The  independent  mode  allows  user  manipulation  of  Centrifuge  Angular  Velocity 
as  a  scalar  from  the  front  panel.  The  default,  integrated  mode  allows  the  input  of  a  one¬ 
dimensional  time  indexed  array  of  Angular  Velocity  linked  from  the  Centrifuge  Dynamic 
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Model.  Visual  feedback  is  employed  to  test  the  logic,  aid  in  troubleshooting,  and  verify 
the  toggle  switch  position. 

2,  Diagram 

a.  Radial  Velocity  of  a  Subject  At-rest 


Logic 

Case  Response 

Function  Start  >  Current  Time 

Function  Duration  +  Function  Start  >  Current  Time 

False 

Function  Start  >  Current  Time 

Function  Duration  +  Function  Start  <=  Current  Time 

True 

Function  Start  <=  Current  Time 

Function  Duration  +  Function  Start  <=  Current  Time 

False 

Table  9:  Motion  Case  Logic 
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The  diagram  of  the  Human  Disorientation  Model  contains  a  Sequence 
Structure  that  consists  of  three  frames.  The  first  frame,  pictured  in  Figure  30,  collects  the 
user’s  inputs  and  computes  radial  velocity  with  respect  to  the  motion  function.  Those 
inputs  are  Radial  Translation  and  Angular  Velocity  scalar  controls  and  the  Human 
Parameters,  Human-Centrifuge  Interaction,  and  Simulation  Parameters  clusters.  The 
motion  function  is  determined  through  the  use  of  a  Ring  Control.  The  default  selection  of 
the  Ring  Control  is  the  Stand  &  Sit  Function.  The  Ring  Control  is  not  designed  to  be 
manipulated  from  a  superior  VI,  therefore  when  the  Centrifuge  Dynamic  Model  is  linked 
within  another  VI,  it  will  generate  its  responses  using  the  Stand  &  Sit  Function  which 
means  that  the  subject  will  stand  and  sit  alternately  in  a  manner  determined  by  the  user. 
In  the  case  of  a  subject  at  rest,  as  shown  in  Figure  30,  the  radial  velocity  will  be  zero  and 
the  radial  position  will  not  change. 

For  each  case,  a  simple  logic  tree,  shown  in  Table  9,  and  Case  Structure  is 
used  to  determine  when  to  apply  the  designated  motion  function.  Radial  Position,  Radial 
Velocity,  and  Time  are  indexed  into  one-dimensional  arrays  and  connected  to  Sequence 
Locals  for  later  use,  bundled  with  one-dimensional  Time  arrays,  and  plotted  versus  Time 
on  the  Front  Panel. 


b.  Radial  Velocity  of  a  Subject  Alternately  Standing  and  Sitting 

The  second  case,  pictured  in  Figure  31,  uses  a  sine  function  to  generate  the 
motion  of  a  person  sitting  and  standing.  The  user  input.  Radial  Translation,  becomes  the 
amplitude  of  the  sine  function.  The  following  equations  were  utilized: 

r  =  (l.55Ar)sin(/At) 


r 


=  r^+  rAt 
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Figure  31:  Determination  of  Subject  Radial  Velocity  when  Alternately  Standing 

and  Sitting 


Figure  32:  Determination  of  Subject  Radial  Velocity  when  Climbing  Up  a  Ladder 
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c.  Radial  Velocity  of  a  Subject  Climbing  Up  a  Ladder 

The  third  case,  pictured  in  Figure  32,  uses  the  negative  absolute  value  of  a 
sine  function  to  generate  the  motion  of  a  person  climbing  up  a  ladder.  The  user  input. 
Radial  Translation,  becomes  the  amplitude  of  that  sine  function.  The  following  sine 
function  was  utilized; 

r  =  -|Arsin(/At)| 


f 


=  r^+  rAt 


d.  Radial  Velocity  of  a  Subject  Climbing  Down  a  Ladder 


Figure  33:  Determination  of  Subject  Radial  Velocity  when  Climbing  Down  a 

Ladder 


The  fourth  case,  pictured  in  Figure  33,  uses  the  absolute  value  of  a  sine 
function  to  generate  the  motion  of  a  person  climbing  down  a  ladder.  The  user  input. 
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Radial  Translation,  becomes  the  amplitude  of  that  sine  function.  The  following  sine 
function  was  utilized; 

r  =  |Arsin(/At)| 


r 


=  r^+  rAt 


e.  Semi-Circular  Canal  Response 

As  with  the  Centrifuge  Steady  State  Response  in  the  Centrifuge  Dynamic 
Model,  the  Semi-Circular  Canal  Response  in  the  Human  Disorientation  Model  is 
computed  within  a  subordinate  VI,  SCC  Hand  Math.  The  various  parameters  of  the 
Semi-Circular  Canal  Model;  Angular  Velocity,  Cupula  Deflection,  Cupula  Velocity, 
Friction  Losses,  Radial  Velocity,  Radial  Position,  and  Restoring  Force  are  linked  to  the 
subordinate  VI  with  the  outputs.  Cupula  Deflection  and  Cupula  Velocity,  linked  to 
Sequence  Locals  for  use  in  a  later  frame. 
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The  use  of  a  subordinate  VI  for  the  computation  of  the  Semi-Circular 
Canal  Response  reduced  program  complexity,  enhanced  troubleshooting,  and  easy 
substitution  of  the  analytical  check  case  for  comparison  purposes. 

f  Display  of  Human  Disorientation  Model  Results 

In  the  final  frame  of  the  Human  Disorientation  Model,  the  results  of  the 
State  Space  Analysis,  Cupula  Deflection  and  Cupula  Velocity,  together  with  Coriolis  and 
G-Level  Experienced  are  bundled  with  a  one-dimensional  Time  array  and  plotted  on  the 
Front  Panel.  Auto-indexing  was  not  used  with  the  controls  in  this  VI  because  when 
running  independently  the  empty  one-dimensional  Angular  Velocity  array  would  conflict 
with  the  auto-indexing  capability  of  LabVIEW. 
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E.  SEMI-CIRCULAR  CANAL  DYNAMIC  MODEL 

1,  Front  Panel 

The  front  panel  of  the  Semi-Circular  Canal  Response,  shown  in  Figure  36,  is  not 
designed  for  direct  user  manipulation.  Rather,  this  VI  is  designed  to  be  a  subordinate  VI 
nested  within  the  Human  Disorientation  State  Space  VI.  The  input  Controls  on  the  Front 
Panel,  which  are  linked  from  the  Human  Disorientation  Model,  populate  the  variables  of 
the  equations  that  compute  the  Semi-Circular  Canal  Response  of  a  person  within  a 
rotating  centrifuge.  The  output  Indicators  on  the  Front  Panel,  Cupula  Deflection  and 
Cupula  Velocity,  link  the  results  of  this  VI  to  the  Human  Disorientation  Model.  The 
major  components  of  this  VI  are  listed  in  Table  10. 


ISCC 

1  Hand 

see  Hand  Math.vi 

Degrees  and  Radians  2  way  Converter.vi 

Table  10:  Major  Components  of  the  Human  Disorientation  Model 


This  VI  Performs  the  Human  Disorientation 

State  Space  Computations 

Sample  Interval,  sec 

y  0.00 

Cupula  Deflection,  deg 

r)  0.00 

Cupula  Deflection,  deg 

Cupula  Velocity,  deg/s 

0.000000 

^  0.00 

Angular  Velocity,  rad/s 
^  0.00 

Restoring  Force  Constant  k/J,  sec 

Cupula  Velocity,  deg/sec 

0.000000 

^  0.00 

Friction  Constant  b/J,  sec 

^  0.00 

Radial  Position,  m 

^  0.00 

Radial  Velocity,  m/s 

^  0.00 

Figure  36:  Semi-Circular  Canal  Dynamic  Model  User  Interface 
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2,  Diagram 

a.  Cupula  Deflection 

The  diagram  of  the  Semi-Circular  Canal  Response  contains  a  Sequence 
Structure  that  consists  of  three  frames.  The  first  frame,  pictured  in  Figure  37,  collects  the 
linked  inputs  and  calculates  Cupula  Deflection.  Cupula  Deflection,  computed  in  radians, 
is  converted  to  and  from  degrees  when  linked  outside  of  this  VI  for  easier  human 
interpretation.  The  third  frame,  not  pictured,  shows  the  system  of  equations  used  in  this 
VI.  The  formula  used  in  this  frame  is  based  on  Equation  (C.45)  in  Appendix  C: 

^2  ;  A  ^  j  1  \  ^2 


^t+M 


1  + 

^  IJ 


V 


-I- 


y 


V 


4J 


4J 


IR 


.  bAt  kAt^ 
1  + - -I-- 


2J  4A 


Figure  37:  Determination  of  Cupula  Deflection 


b.  Cupula  Velocity 

The  second  frame,  pictured  in  Figure  38,  collects  the  linked  inputs  and 
calculates  Cupula  Velocity.  Cupula  Velocity,  computed  in  radians  per  second,  is 
converted  to  and  from  degrees/second  when  linked  outside  of  this  VI  for  easier  human 
interpretation.  The  formula  used  in  this  frame  is  based  on  Equation  (C.46)  in  Appendix 
C: 


49 


^  k^t^( ^  (pAt 


^i+isi 


IJ 


V 


^  ^  kip^At  ■  IcoVj^At'^ 
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Compute  Cupula  Velocity  (phi  dot)  based 
on  current  state  and  predict  new  state 


numi  =(kJ’‘dt)/'2; 
num2=phi+((phidot’‘dl)/2); 
num3=(kJ’‘phi’‘dt)/2; 
num4=(bJ  ’■phidot’'dt)/2; 
num5=(2’'omega’'v’'dl)/fj 
den1=(bJ’'dt)/2; 
den2=(U1dr-2))/'4; 

:oniega|  phidon=((-num1  ’‘num2)+(-num3+phidot-num4+num5))/' 
(1  +den1  +den2); 
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Figure  38:  Determination  of  Cupula  Velocity 
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IV.  RESULTS 


A,  7.5  FOOT  RADIUS  CENTRIFUGE  FOR  VALIDATION  EXPERIMENT 


Figure  39:  Experimental  Apparatus  and  Volunteer  Test  Subjects 


A  7ft.  6.5in.  (1.15m)  diameter  Merry-Go-Round,  pietured  in  Figure  39  was  found 
on  the  Former  Fort  Ord  in  Marina,  CA  that  enabled  the  eonduet  of  a  non-seientific 
validation  experiment.  This  brief  experiment,  eondueted  between  9  and  10  AM  on  10 
September  2002,  served  to  determine  the  relationship  between  the  computer  model’s 
quantitative  results  and  a  person’s  qualitative  experience.  This  experiment  should  be 
viewed  as  a  rough  guide  for  those  conducting  follow-on  studies  to  establish  the 
parameters  to  scientific  experiments  or  the  analysis  of  preexisting  data. 

1,  Test  Subjects 

A  small  number  of  test  subjects  were  selected  from  Naval  Postgraduate  School 
students  in  the  Space  Systems  Operations  curriculum  graduating  in  September  2002. 
They  were  selected  due  to  their  acquaintance  with  the  author,  their  scheduling  flexibility, 

and  the  variety  of  their  military  experience.  The  volunteers  were  briefed  on  how  the 
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experiment  would  be  eondueted,  and  were  instrueted  to  report  the  diseomfort  they 
experieneed  during  the  three  phases  of  eaeh  run.  The  small  population  of  test  subjeets 
used  in  this  experiment  do  not  eonstitute  a  seientifie  sampling  of  results;  however,  these 
results  provided  aneedotal  evidenee  of  the  level  of  diseomfort  experieneed  at  partieular 
eupula  veloeities. 


Name 

Baekground 

LT  Chris  Howse,  USNR 

Naval  Surfaee  Warfare,  Destroyers 

CPT  Mike  York,  USA 

Army  Air  Defense  Artillery,  Patriot 

LT  Wes  Sanders,  USN 

Naval  Flight  Offieer,  E-2C  Hawkeye 

Table  11:  Test  Subjects  and  Their  Backgrounds 


2,  0  Revolutions  Per  Minute 

Cyelieal  movement  of  the  head  in  a  normal,  non-rotating,  environment  did  not 
produee  disorientation  in  any  of  the  test  subjeets. 

3.  10  Revolutions  Per  Minute 

a.  Model  Predictions 

The  LabView  eomputer  model  developed  using  the  analysis  eontained  in 
Appendix  C  predieted  Cupula  Defleetions  ranging  from  -5°  to  1.75°  and  Cupula 
Veloeities  as  great  as  ±9°/see.  The  predieted  Centrifugal  Foree  was  a  maximum  of  about 
0.13G. 

b.  Subject  Experiences 

Utilizing  the  disorientation  seale  developed  in  Table  3,  the  test  subjeets 
were  exposed  (as  elosely  as  available  equipment  would  permit)  to  the  modeled 
eonditions.  During  the  first  phase  all  subjeets  reported  no  feeling  of  disorientation  or  a 
Disorientation  Level  of  1  (Level  1).  In  the  seeond  phase  all  subjeets  reported  a  Level  2. 
In  the  third  phase  all  subjeets  reported  Level  1  to  1.5. 
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This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 


see  eharacteristics 
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Figure  40:  10  RPM  Validation  Experiment  Prediction 


4,  20  Revolutions  Per  Minute 

a.  Model  Predictions 

The  LabView  computer  model  predicted  Cupula  Deflections  ranging  from 
-10°  to  3°,  Cupula  Velocities  as  great  as  ±18°/sec,  and  Centrifugal  Force  at  a  maximum 
of  about  0.51  G. 

b.  Subject  Experiences 

Once  again  using  the  scale  developed  in  Table  3,  the  test  subjects  were 
exposed  to  the  modeled  conditions.  During  the  first  phase  all  subjects  reported  a  Level  2 
sensation  largely  resulting  from  Centrifugal  Force.  In  the  second  phase,  a  Level  4  feeling 
was  reported  with  the  addition  of  an  oculogyral  illusion  of  the  horizon  rolling  in  phase 
with  the  head  motion,  a  direct  indication  of  the  Coriolis  induced  torque  in  the  axial  plane 
shown  in  Equation  (C.IO).  Finally,  in  the  third  phase  a  Level  3  feeling  was  reported. 
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Figure  41:  LT  Sanders  Demonstrating  the  Test  Subject’s  Perspective 


Figure  42:  LT  Sanders  Demonstrating  Experiment  Phases  One  and  Three, 


Figure  43:  The  Author  Demonstrating  Maximum  Radial  Position  of  the  Head 
During  Phase  Two,  Photograph  by  LT  Wesley  Sanders,  USN 
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Figure  44:  The  Author  Demonstrating  Minimum  Radial  Position  of  the  Head, 
Photograph  hy  LT  Wesley  Sanders,  USN 

This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 
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Figure  45:  20  RPM  Validation  Experiment  Prediction 


c.  Discussion  of  Validation  Experiment  Results 
Of  particular  note  is  the  direct  relationship  between  Angular  Veloeity, 
Coriolis,  and  Cupula  Veloeity  shown  in  the  computer  model.  This  relationship  was 
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repeated  during  the  validation  experiment  as  shown  below  in  Table  12.  It  is  also  worth 
noting  that  Centrifugal  Force  quadrupled  when  Angular  Velocity  was  doubled. 


10  RPM 

20  RPM 

Phase  1 

1 

2 

Phase  2 

2 

4 

Phase  3 

1  to  1.5 

3 

Cupula  Velocity 

±9°/sec 

±18°/sec 

Centrifugal  G 

0.13 

0.52 

Table  12:  Validation  Experiment  Results 


Using  the  results  in  Table  12,  larger  centrifuges  can  be  modeled  by 
changing  the  parameters  of  the  centrifuge.  Noting  that  a  minimal  amount  of  discomfort 
was  experienced  at  10  RPM  on  a  1.15  m  Centrifuge.  The  computer  model  shall  evaluate 
1.25  m  (ISS  CAM),  2.5  m,  5  m,  and  10  m  Centrifuges  at  no  more  than  ±9°/sec  of  Cupula 
Velocity. 


B,  1,25  METER  RADIUS  CENTRIFUGE 

The  ISS  Centrifuge  (not  much  larger  than  the  Merry-Go-Round  used  in  the 
validation  experiment)  is  designed  to  accommodate  small  biological  experiments  within  a 
gravitational  environment  that  ranges  from  0.01  G  to  2  G^®.  The  ISS  Centrifuge  is  not 
designed  to  accommodate  humans.  However,  should  experiments  be  carried  out  to  study 
the  ability  of  a  subject  to  tolerate  Coriolis  induced  disorientation,  the  results  in  Figure  46 
should  serve  to  guide  their  development.  The  computer  model,  as  shown  in  Figure  46, 
predicts  Cupula  Deflection  results  equivalent  to  those  in  the  lORPM  Validation 
Experiment  will  be  achieved  at  11  RPM  and  0.17  G  when  the  same  head  motion  is 
performed  at  the  centrifuge’s  maximum  radius.  Note  that  10  cm  of  increased  centrifuge 
radius  yields  0.04  G  of  increased  Centrifugal  Force  at  a  level  of  disorientation  that  was 
deemed  tolerable  by  the  test  subjects  in  the  Validation  Experiment. 

“Centrifuge” 
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This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 

Radial  VelodV,  m/sec 
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Figure  46:  1,25  m  Centrifuge  Human  Disorientation  Prediction 


Additionally,  the  ISS  Centrifuge  Speeifieations  shown  in  Table  1  allow  for  the 
development  of  a  Centrifuge  Dynamic  Model  for  the  purposes  of  scaling  the  properties  of 
the  ISS  Centrifuge  to  fit  larger  centrifuges  more  suited  for  human  occupation.  Figure  47 
details  the  results  of  the  Centrifuge  Dynamic  Model  when  the  ISS  Centrifuge 
Specifications  are  used  as  the  input  parameters  (2  G,  1 .25  m  throw  arm/inner  radius,  and 
300  sec.  spin-up). 

Increasing  the  EMF  Proportionality  Constant  so  that  the  Centrifugal  Force  is  1  G 
yields  the  results  shown  in  Figure  48.  By  leaving  the  Electrical  Parameters  constant  and 
changing  the  Inner  and  Outer  Radius,  the  Centrifuge  Dynamic  Model  can  be  scaled  to 
approximate  the  effects  of  a  similar,  but  larger,  centrifuge  on  a  spacecraft  like  the  ISS. 
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Number  Of  Points 
40001 

This  VI  computes  the  dynamic  response 
of  a  centrifuge  based  on  an  arbitrary  input  Current  Time  Point 
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Figure  47:  Centrifuge  Dynamic  Model  of  the  ISS  Centrifuge  at  2  G 
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This  VI  computes  the  dynamic  response 
of  a  centrifuge  based  on  an  arbitrary  input  Current  Time  Point 
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Figure  48:  Centrifuge  Dynamic  Model  Of  The  ISS  Centrifuge  At  1  G 
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C.  2.5  METER  RADIUS  CENTRIFUGE 

A  centrifuge  with  a  2.5  meter  radius  eould  contain  a  compartment  large  enough  to 
aeeommodate  a  person  and  some  limited  exereise  equipment.  What  would  be  the  effeet 
of  Coriolis  on  a  person  exereising  within  sueh  a  rotating  environment.  Before  modeling 
these  eonditions,  the  following  assumptions  must  first  be  made:  the  astronaut  is  1.83  m. 
(6  ft.)  tall,  when  moving  from  a  seated  to  a  standing  position  the  astronaut  will  translate 
radially  through  0.3  m  (1  ft.)  to  simulate  a  normal  range  of  movement.  The  astronaut  will 
start  from  a  squatting  position  with  his  head  at  a  radial  position  of  0.97  m  (as  if 
performing  military  presses  on  an  exereise  maehine),  the  astronaut’s  feet  will  be  at  a 
radial  position  of  2.5  m,  and  the  eentrifuge  will  be  spinning  at  18.9  RPM  to  attain  1  G  at 
the  iimer  radius  of  the  rotating  drum. 


This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 
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Figure  49:  2,5  Meter  Centrifuge  Human  Disorientation  Prediction  at  1  G 


A  Human  Disorientation  Model  of  the  2.5  meter  centrifuge  at  1  G  predieted  a 
Cupula  Veloeity  of  ±70°/sec  and  Cupula  Deflections  between  -20°  and  10°.  These 
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results  are  significantly  greater  than  the  Cupula  Velocity  and  Defiection  experienced 
during  the  20  RPM  Validation  Experiment  and  would  likely  make  the  astronaut 
experiencing  them  very  ill.  The  results  are  pictured  in  Figure  49. 

Reducing  the  centrifuge  angular  velocity  to  2.5  RPM  yields  results  that  more 
closely  match  the  results  from  the  validation  experiment.  As  shown  in  Figure  50,  Cupula 
Velocity  was  ±9°/sec  and  Cupula  Deflection  ranged  between  1.75°  to  -5°.  However,  the 
Centrifugal  Force  at  centrifuge’s  inner  radius  was  less  than  0.02  G.  Based  on  these 
results,  a  2.5  Meter  Centrifuge  would  not  provide  a  suitable  environment  for  an  astronaut 
to  move  about  normally  and  exercise. 


This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 
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Figure  50:  2,5  Meter  Centrifuge  with  Minimal  Disorientation 


D.  5  METER  RADIUS  CENTRIFUGE 
1.  Human  Disorientation 

What  would  be  the  effect  of  Coriolis  on  a  person  moving  about  within  a  5  meter 


radius  centrifuge?  Perhaps  such  a  centrifuge  could  contain  an  exercise  facility  for 
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astronauts  on  a  long  term  mission.  Adapting  the  assumptions  from  the  2.5  m  Centrifuge 
Analysis,  the  astronaut  is  1.83  m  (6  ft)  tall.  When  moving  from  a  seated  to  a  standing 
position  the  astronaut  will  translate  through  0.3  m  in  the  radial  plane,  start  from  a 
squatting  position  with  his  head  at  a  radial  position  of  3.47  m,  and  his  feet  will  be  at  a 
radial  position  of  5  m.  The  centrifuge  will  be  spinning  at  13.38  RPM  to  attain  1  G  at  the 
inner  radius  of  the  rotating  drum. 

In  the  5  m  centrifuge.  Cupula  Velocity  at  1  G  was  predicted  to  be  ±12°/sec  and 
the  Cupula  Deflection  results  ranged  between  -7°  and  2.25°.  These  results  are  within  the 
range  of  tolerance  experienced  during  the  Validation  Experiment  and  an  astronaut  could 
possibly  adapt  to  their  effects.  The  results  are  pictured  in  Figure  5 1 . 


This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 
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Figure  51:  5  Meter  Centrifuge  Human  Disorientation  Prediction  at  1  G 


Reducing  the  centrifuge  angular  velocity  to  10  RPM  yields  results  that  more 
closely  match  the  results  from  the  10  RPM  validation  experiment.  As  shown  in  Figure 
52,  Cupula  Velocity  was  ±9°/sec  and  Cupula  Deflection  ranged  between  -5°  to  1.5°. 
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However,  the  G-Level  at  the  outer  radius  was  reduced  to  0.56  G.  Based  on  these  results, 
a  5  Meter  Centrifuge  would  provide  an  astronaut  the  benefit  of  “Fractional  G”  to  move 
about  and  exercise.  However,  at  1  G,  there  would  be  some  disorientation  and  its 
usefulness  would  depend  on  the  astronaut’s  ability  to  adapt. 

2,  Centrifuge  Dynamic  Model 

Because  of  the  possibility  that  a  5  m  Centrifuge  producing  1  G  of  Centrifugal 
Force  might  be  useful;  the  Centrifuge  Dynamic  Model  in  Figure  48  was  reconfigured  to 
have  an  Outer  Radius  of  5.01  m,  an  Inner  Radius  of  5.00  m,  and  an  EMF  Proportionality 
Constant  of  19.9.  The  results  pictured  in  Figure  53  show  that  such  a  Centrifuge  would 
require  over  1  hour  to  spin-up  and  would  weigh  about  20200  kg. 


This  VI  computes  the  Semi-Circular  Canal  response  for  an  arbitrary  stimulus  inside  a  centrifuge 
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Figure  52:  5  Meter  Centrifuge  with  Minimal  Disorientation 
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Number  Of  Points 
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This  Vi  computes  the  dynamic  response 
of  a  centrifuge  based  on  an  arbitrary  input  Current  Time  Point 
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Figure  53:  Centrifuge  Dynamic  Model  of  a  5  Meter  Centrifuge  at  1  G 


E.  10  METER  RADIUS  CENTRIFUGE 

1.  Human  Disorientation 

Again,  what  would  be  the  effect  of  Coriolis  on  a  person  moving  within  a  1 0  meter 
radius  centrifuge?  Perhaps  such  a  centrifuge  could  serve  as  a  complete  habitat  containing 
berthing,  exercise,  and  food  preparation  facilities.  Adapting  the  assumptions  from  the  2.5 
m  Centrifuge  Analysis,  the  astronaut  is  1.83  m  (6  ft)  tall.  When  moving  from  a  seated  to 
a  standing  position  the  astronaut  will  translate  through  0.3  m  in  the  radial  plane,  start 
from  a  squatting  position  with  his  head  at  a  radial  position  of  8.47  m,  and  his  feet  will  be 
at  a  radial  position  of  10  m.  The  centrifuge  will  be  spinning  at  9.46  RPM  to  attain  1  G  at 
the  inner  radius  of  the  rotating  drum. 

In  the  10  meter  centrifuge  at  1  G,  Cupula  Velocity  was  predicted  to  be  ±2.75°/sec 
and  Cupula  Deflection  ranged  between  -1.6°  and  0.6°.  These  results,  pictured  in  Figure 
54,  are  less  than  those  from  the  10  RPM  validation  experiment  and  should  be  well 
tolerated  by  an  astronaut.  Based  on  these  results,  a  10  Meter  Centrifuge  would  provide 
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an  astronaut  the  benefit  of  a  1  G  gravity  field  to  move  about  and  exereise  without 
signifioant  disorientation. 


This  VI  computes  the  Semi-Circu 

see  Characteristics 

Human-Centrifuge  Interaction 

Restoring  Force,  k/J,  I/sec 

Initial  Radial  Position,  m 

rj!20.00 

^8.47 

Friction  Losses,  b/J,  Vsec 

Function  Speed,  Hz 

;j|2oaoo 

:j(a50 

Simulation  Parameters 

Human  Parameters 

Input  Start  Time,  sec 

Function  Type 

r )  Stand  &  Sit  1 

Input  Start  Function  Duration,  sec  Subi.  Function  Start  Time,  sec 

^/aoo 

yi  60.00 

Sim  Stop  Time,  sec 

Subiect  Function  Duration,  sec 

:^lfi20.od~ 

:jfio.oo 

Sim  Sample  Interval,  sec 

Radial  Translation,  m 

0,01 000 

;_)  0.3000 

Input  Stop  Time,  sec 

Angular  Velocity,  RPM 

;)  9.4800 

12000 

Number  Of  Points 
12001 

Current  Time,  sec 
120.00 

Angular  VelociV,  rad/sec 
0.9906 

Independent 

i 

Integrated 


Coriolis,  N 
100.0- 


HHIHinilHI 

HIIIIIIIIIIHI 

1.000  55.000  60.000  65.000  70.000  75.000  80.000  85.000  90.C 

G-Level  Experienced,  G 
0.850-1 
„  0.840- 
1  0.830- 
£  0.820- 


0.810 

50.000  55.000 


HUnirWIIIHHHH 


60.000  65.000  70.000  75.000  80.000  85.000  90.000 
Time 


Radial  VelociV-  m/sec 
50.0-1 
25.0-1 

u 

I  0-0- 

I  -25.0- 


-50.0-' 


jinifiMii 

ininniin 

iniiiiiiiii 

iwfliwr 


50.0  55.0  60.0  65.0 


70.0 

Time 


75.0  80.0  85.0  90.0 


Radial  Position,  m 


1.5- 

1.4- 

HIIIIIIIIIIIIH 

HnniiiiHH 

1.0  55 

'.0  60.0  650  70.0  750  8C 

i.O  85 

i.O  90. 

Time 

Cupula  VelodV-  Vsec 


HiinuimH 

HI1IIIII1IH 

lOOO  55000  60,000  65.000  70.000  751 

300  80.1 

300  85.1 

300  90.C 

Cupula  Deflection,  * 

1.0-H 


•2.0-j 

50.000 


IlMMlilll 

— vmmir 


55.000  60.000  65.000 


70.000 

Time 


75.000  80.000  85.000  90.000 


Figure  54:  10  Meter  Centrifuge  Human  Disorientation  Prediction  at  1  G 


2,  Centrifuge  Dynamic  Model 

Because  of  the  possibility  that  a  5  m  Centrifuge  producing  1  G  of  Centrifugal 
Force  might  be  useful;  the  Centrifuge  Dynamic  Model  in  Figure  48  was  reconfigured  to 
have  an  Outer  Radius  of  10.01  m,  an  Inner  Radius  of  10.00  m,  and  an  EMF 
Proportionality  Constant  of  28.2.  The  results  pictured  in  Figure  53  show  that  such  a 
Centrifuge  would  require  nearly  19.5  hours  to  spin-up  and  would  weigh  about  70800  kg. 


F.  COLLECTED  RESULTS 

The  results  of  the  model  centrifuges,  with  respect  to  the  human  tolerance  of 
Coriolis  found  during  the  Validation  Experiment  are  listed  in  Table  13. 
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This  VI  computes  the  dynamic  response 
of  a  centrifuge  based  on  an  arbitrary  input 
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^28.200 

3849193. 

^  Step  0 

Centrifuge  Dimensions 

Simulation  Parameters 

Outer  Radius,  m 

Input  Start  Time,  sec 

:j|o.oo 

Inner  Radius,  m 

;)lio.oo 

Cap  Thickness,  m 

Shaft  Mass  Density,  kg/nf 
;j|l0000.0C 

Input  Start  Function  Duration,  sec 

;jro.oo 

yjo.01 

Shaft  Length,  m 
'J|3.00 

Shaft  Radius,  m 

Centrifuge  Mass  Density,  kg/nf 
:j|l0000.0C 

Sim  Stop  Time,  sec 

;_^'8o6ooir 

Sim  Sample  Interval,  sec 

^^(oToooo^ 

Drum  Length,  m 

Input  Stop  Time,  sec 

:^80000.0C 

Figure  55:  Centrifuge  Dynamic  Model  of  a  10  Meter  Centrifuge  at  1  G 


Radius 

m 

Radial 

Position 

m 

Tolerable 

Speed 

RPM 

Tolerable 

Centrifugal 

Force 

RPM 

Spin  Up 
Time 
min 

Mass 

kg 

Peak 

Power 

W 

Steady  State 
Current 

A 

Power 

Required 

mW/s 

2.5 

0.97 

2.5 

0.02 

1.15 

1.15 

10 

0.13 

1.25 

1.25 

11 

0.17 

5 

2807 

150 

0.047 

43.7 

5 

3.47 

10 

0.54 

133 

20173 

150 

0.023 

1.936 

10 

8.47 

9.46 

1 

1167 

70815 

150 

0.016 

0.282625 

Table  13:  Collected  Results 


The  relationship  between  the  radial  position  of  the  head  and  tolerable  RPM  is 
shown  in  Figure  56.  Each  point  is  labeled  with  the  radial  position  associated  with  it.  The 
relationship  between  the  radial  position  of  the  head  and  tolerable  RPM  is  shown  in  Figure 
57.  Each  point  is  labeled  with  the  radial  position  associated  with  it.  The  relationship 
between  the  centrifuge  radius  and  spin-up  time  is  shown  in  Figure  58.  Each  point  is 
labeled  with  the  centrifuge  radius  associated  with  it.  The  relationship  between  centrifuge 
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radius  and  centrifuge  mass  is  shown  in  Figure  59.  Each  point  is  labeled  with  the 
centrifuge  radius  associated  with  it. 


Figure  56:  Tolerable  Angular  Velocity  Versus  Radial  Position 


Centrifugal  Force 


Figure  57:  Tolerable  Centrifugal  Force  Versus  Radial  Position 
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Figure  59:  Centrifuge  Mass  Versus  Centrifuge  Radius 
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V.  CONCLUSION 


A.  VALIDATION  EXPERIMENT 

1.  Effects  of  Coriolis  on  Human  Disorientation 

Movement  of  the  head  in  the  radial  plane  of  a  rotating  environment  generates  a 
profound  feeling  of  disorientation,  an  oeulogyral  illusion  of  the  horizon  twisting,  and 
depending  on  the  angular  veloeity  eould  induee  motion  siekness.  As  demonstrated 
during  the  validation  experiment,  these  feelings  are  direetly  proportional  to  the  angular 
veloeity  of  the  eentrifuge.  The  validation  experiment  showed  that  slower  angular 
veloeities  generate  less  disorientation. 

2.  Validation  Experiment 

The  validation  experiment  was  designed  to  determine  how  the  physieal  response 
of  the  Semi-Cireular  Canals,  as  reprodueed  by  the  Human  Disorientation  Model, 
eorresponds  to  the  various  feelings  assoeiated  with  disorientation  in  a  rotating 
environment.  Utilizing  a  Merry-Go-Round  to  provide  a  1.15  m  Rotating  Platform,  three 
subjeets  were  subjeeted  to  rotational  motion  at  10  and  20  RPM.  At  10  RPM,  a  very  mild 
sense  of  disorientation,  akin  to  slight  dizziness,  was  experieneed  when  eaeh  subjeet 
moved  their  head  baek  and  forth,  radially,  through  10  em  for  10  see.  No  disorientation 
was  experieneed  at  10  RPM  during  spin  up,  spin  down,  and  steady  state  rotation  without 
head  movement.  At  20  RPM,  a  strong  sense  of  disorientation,  akin  to  nausea,  was 
experieneed  aeeompanied  by  an  illusion  of  the  horizon  tilting  as  eaeh  subjeet  moved  their 
head  baek  and  forth,  radially,  through  10  em  for  10  see.  The  illusion  moved  in  phase 
with  the  head  movement.  A  mild  sense  of  disorientation  was  experieneed  at  20  RPM 
during  spin  up,  spin  down,  and  steady  state  rotation  that  eaeh  subjeet  related  to  their 
efforts  to  maintain  their  position  against  the  strong,  0.52  G  Centrifugal  Foree. 

3.  Human  Tolerance  of  Coriolis 

The  Human  Disorientation  Model  presented  in  this  paper  shows  that  the  ability  of 
a  person  to  tolerate  Coriolis  is  direetly  related  to  Cupula  Veloeity  and  not  Cupula 
Defleotion.  Cupula  Defleetion  other  than  zero  stimulates  the  vestibular  nerve  the 
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magnitude  of  which  is  directly  interpreted  by  the  brain  as  rotational  motion  in  the  plane 
of  the  associated  Semi-Circular  Canal.  Cupula  Velocity  is  dependent  on  the  torque 
produced  during  radial  translation.  This  torque  is  the  vector  product  of  Coriolis  Force 
and  Radial  Velocity.  The  Human  Disorientation  Model  showed  that  Cupula  Velocity  for 
the  same  Radial  Motion  and  Centrifugal  Force  was  reduced  when  Centrifuge  Radius  was 
increased  and  Angular  Velocity  was  decreased. 

Assuming  a  1  G  Centrifugal  Force,  a  centrifuge  with  a  5  m  radius  is  the  smallest 
that  could  be  occupied  by  an  astronaut  where  the  astronaut  could  withstand  the  effects  of 
disorientation.  However,  a  5  m  centrifuge  imparts  enough  disorientation  that  motion 
sickness  would  be  experienced,  impairing  the  mission.  The  problems  associated  with 
disorientation  do  not  exist  in  the  10  m  centrifuge. 

4,  Man  Rated  Centrifuge  Characteristics  Given  Peak  Electrical  Power 
Load 

This  report  calculated  the  dynamic  characteristics  of  5  and  10  m  centrifuges  based 
on  a  150W  Peak  Load  primer  mover.  The  5  and  10  m  centrifuges  simulated  with  the 
human  dynamic  model  were  input  into  the  Centrifuge  Dynamic  Model  of  the  ISS 
Centrifuge.  Fixing  the  characteristics  of  Steady  State  Centrifugal  Force  and  Peak  Start- 
Up  Power  required  only  the  EMF  Proportionality  Constant  be  changed  together  with  the 
Inner  and  Outer  Radius.  The  results  were  that  a  10  m  radius  centrifuge  would  spin  up  in 
more  than  19  hours  and  weigh  70815  kg,  making  it  unfeasible  in  the  near  term.  The  5  m 
radius  centrifuge,  on  the  other  hand,  weighing  in  at  20173  kg  with  a  133  min  spin-up  time 
is  a  feasible  project  that  merits  further  consideration. 


70 


VI.  FOLLOW-ON  RESEARCH 


A.  GROUND-BASED  RESEARCH  OPPORTUNITIES 

1.  Correlation  with  Existing  Data 

Existing  data  from  previous  human  centrifuge  studies  can  be  utilized  to  further 
evaluate  the  model  presented  in  this  paper.  Additionally,  scientific  research  can  be 
designed  to  further  evaluate  disorientation  in  a  manner  that  best  simulates  the  space 
environment. 

2.  Manned  Centrifuge  Research 

Existing  centrifuges  in  the  military  and  civilian  sector  can  be  utilized,  together 
with  these  models,  to  conduct  scientific  research  into  the  level  of  disorientation 
experienced  compared  to  the  physical  response  of  the  Semi-Circular  Canals  under 
differing  conditions  in  a  rotating  environment. 

3.  Dual-Use  of  a  Manned  Centrifuge  as  a  Momentum  Wheel 

What  are  the  benefits  of  a  centrifuge  when  it  is  also  used  as  a  momentum  wheel 
for  the  purpose  of  maintaining  spacecraft  stability?  Would  it  be  beneficial  to  use  a 
centrifuge  to  replace  control  moment  gyroscopes  or  standard  momentum  wheels? 

4.  Dynamic  Properties  of  Other  Centrifuge  Types  Versus  Disorientation 

Given  that  Human  Disorientation  in  a  rotating  environment  is  independent  of  the 
design  of  the  centrifuge,  would  a  radial-arm  or  tethered  centrifuge  be  more  feasible  for 
deployment  in  the  near-term?  Could  one  be  installed  as  an  upgrade  to  the  ISS  or  a  future 
space  station? 

B,  POTENTIAL  SPACE-BASED  RESEARCH 

1.  Use  of  The  Disorientation  Model  to  Design  ISS  Centrifuge  Research 

How  will  the  animals  proposed  for  use  in  ISS  Centrifuge  research  respond  to  the 
rotating  environment  within  that  centrifuge?  The  Human  Disorientation  model  allows  for 
the  Eriction  Constant  and  Restoring  Eorce  Constant  of  the  cupula  within  any  Semi- 
Circular  Canal  to  be  adjusted  to  fit  those  of  any  animal  with  vestibular  organs  similar  to 
those  in  humans. 
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VII.  APPENDICES 


A,  APPENDIX  A 


1,  Modeling  the  External  Dynamics  of  a  Centrifuge 

This  section  discusses  the  mathematical  rigor  behind  the  dynamics  of  a  centrifuge. 
Equation  Chapter  (Next)  Section  1 


Analytical  Term 

Meaning 

Units 

B 

Friction  Fosses  due  to  supports,  bearings,  and  gears 

N 

i 

Motor  Current 

A 

J 

Centrifuge  Moment  of  Inertia 

kg  m^ 

kf 

EMF  Coefficient 

Unitless 

K 

Torque  Coefficient 

Unitless 

L 

Motor  Inductance 

H 

R 

Motor  Resistance 

kl 

V 

Drive  Voltage 

V 

Foad  Induced  Back  EMF 

V 

CO 

Centrifuge  Angular  Velocity 

rad/s 

d 

Centrifuge  Angular  Position 

rad 

T 

Motor  Torque 

N-m 

Table  14:  Centrifuge  External  Terms 


2,  Analytical  Model 

The  external  analysis  of  the  effects  of  a  centrifuge  (as  illustrated  in  Figure  60) 
upon  the  spacecraft  where  it  is  installed  is  important  in  determining  the  size  of  centrifuge 
that  optimizes  human  comfort  and  spacecraft  bus  impact. 
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Figure  60:  Centrifuge  External  Terms 


a.  Transfer  Functions 

The  differential  equations  deseribing  the  instantaneous  eurrent  flow  within 
the  Prime  Mover  and  the  Centrifuge  are^^; 


di  ,  do 
Ri  +  L  —  =  v-kf  — 
dt  dt 

,  .  ^d^O  ^dO 
kd=J — ^  +  B — 
dt  dt 


(A.1) 

(A.2) 


The  Laplaee  Transform  of  equations  (A.l)  and  (A.2)  are: 

{R  +  sL)I  =  V-kfS@  (A.3) 

kJ  =  {Js^  +Bs)@  (A.4) 

Solving  both  equations  with  respeet  to  eurrent  in  the  Laplaee  Domain: 


(A.5) 


{R  +  sL)  k^ 

Solving  for  angular  veloeity  returns  the  transfer  funetion  that  relates 
angular  veloeity  to  the  drive  voltage: 


John  J.  D’Azzo  and  Constantine  H.  Houpis,  Feedback  Control  System  Analysis  and  Synthesis.  2"'^ 
Ed.,  New  York:  McGraw  Fiill,  1966,  pp.  44-47 
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co{s)  =  s&  = 


-V 


(A.6) 


[^(i?  +  sL)  [Js  +  B)  +  j 

Substituting  equation  (A.6)  back  into  equation  (A.5)  returns  the  transfer 
function  that  relates  motor  armature  current  to  the  drive  voltage: 

/  (Js  +  B) 

-  =  T - - - - n  (A.7) 

V  [{R  +  sL){Js  +  B)  +  k^kf] 

The  instantaneous  power  consumed  by  the  centrifuge  and  its  prime  mover 


is  described  by: 


p  =  I.V  =  V" 


{Js  +  B) 


^{R  +  sL){Js  +  B)  +  k^k 

b.  Steady  State  Response 
Letting: 

cd{s)  =  s& 

Equation  (A.6)  can  be  rewritten  as: 


(A.8) 


(A.9) 


0)(s)  ■ 


(kA 

{lj) 


2  , 
5  + 


BL  +  RJ 
LJ 


5  + 


f  BR  +  k^kA 

LJ 


K(s) 


(A.10) 


If  the  input  voltage  is  assumed  to  be  a  step  input,  then  the  Laplace 
Transform  of  that  input  voltage  is  (the  subscript  ‘ss’  represents  steady  state): 

V{s)  =  ^ 

s 

Substituting  for  the  input  voltage  yields: 


(A.11) 


0){s)  ■■ 


(kA 

{lj) 


2  , 
5  + 


BL  +  RJ 
LJ 


5  + 


f  BR  +  k^kA 

LJ 


V  s  J 


(A.12) 


velocity: 


Utilizing  the  final  value  theorem  to  obtain  the  steady  state  angular 


a  =  limf^tyf^)) 
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(A.13) 


T7^ - u\  r  r 

[R  +  sL)[Js  +  B)  +  k^k^ 


{RB+\k/^^ 


(A.14) 


(A.15) 


Substituting  the  results  from  equation  (A.15)  into  equation  (A.5)  yields  the 
steady  state  eurrent  required  to  maintain  a  eonstant  angular  veloeity: 


hs  r  ^SS 

K 


(A.16) 


(RB  +  k,k^)' 

Steady  state  power  consumption  is: 


P  =  V  i 

ss  ss  ss 


(A.17) 


(A.  18) 


P  = - =- - 

”  (RB  +  k,k^) 


(RB  +  k,kf) 


(RB  +  k,k^f 


(A.19) 


(A.20) 


(A.21) 


c.  Transient  Response 

Determining  the  centrifuge’s  response  to  transient  events  is  necessary  to 
determine  it’s  Engineering  Cost  and  effect  on  its  human  occupants.  Once  again, 
assuming  a  step  input  for  the  input  voltage,  the  angular  velocity  transfer  function  can  be 


written  as: 


(BR  +  k^k^) 


Ks  +  ^^2  fBL  +  Rj]^  ^  r BR  +  k,kf 


(A.22) 
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The  dynamic  response  is  contingent  upon  the  roots  of  the  following 
characteristic  equation: 


2  , 
5  + 


BL  +  RJ 


5  + 


BR  +  k^k 


T-'-f 


LJ 


=  0 


(A.23) 


V  LJ 

^  'V  / 

The  characteristic  equation  in  equation  (A.23)  is  in  a  form  that  shows  the 
system  has  one  or  more  complex  conjugate  roots.  Therefore  the  damping  ratio  (^)  and 
natural  frequency  ((Dn)  can  be  determined  using  the  following  equation's ; 

s^  +  2,^0)„s  +  0)l 

f  BL  +  RJ 


LJ 


,0)^ 


BR  +  kJif^ 


V 


LJ 


(A.24) 

(A.25) 


The  complex  roots  of  equation  (A.24)  are: 


(A.26) 


Defining  a  ‘damped  natural  frequency  ((Dd)’  as: 

The  roots  in  equation  (A.26)  may  be  simplified  as: 

sx=-^a>n^a>d 

S2=-^0)n-0), 

Using  equation  (A.28),  equation  (A.22)  may  be  rewritten  in  terms  of 
damping  ratio  and  natural  frequency  as  follows: 

(A.29) 


(A.27) 


(A.28) 


K  ] 

(  2  ^ 

OJn 

1  V^s  J 

\\BR  +  Kkf)^ 

^s(s^+2<^0)„s  +  0)l)^ 

^  co{s)^ 

= 

K  ] 

f  (BR  +  k,k^)"\ 

1 

1  Ks  J 

{{BR  +  k,kf)^ 

JL 

V  y 

^60(s)^ 

1 

V  Ks  y 

UlJ 

[5(5-5, )(5-52)J 

(A.30) 


(A.31) 
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D’Azzo  and  Houpis,  p.  246  and  equation  7-110 
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(A.32) 


Utilizing  Laplace  Transform  Tables^^^  the  inverse  Laplaee  Transform  of 
equation  (A.32)  is: 


=  ^  £-‘  - 
[JLJ  s 


(s  +  (^co„-co,}}(s  +  (^co„-co,}} 


F(s)  = 


(s  +  0}(s  +  (^co„-co,}}(s  +  (^co„-co,}} 


(A.33) 


(A.34) 


—ai  —ui  —Cl 

{b-a)[c-a)  [c-b)[a-b)  [a-c){b-c) 


m= 


{<^CO„+CO,)-{<^CO„-CO,)  J  [(^£0„+£0,)-(^60„-£0,) 


Equation  (A. 3 6)  reduees  to: 


(i«. -(i«. 

ft».v  2 


By  noting  that: 


(A.35) 


(A.36) 


(A.37) 


)'  -  OJ]  ^^0)1  -  (ol  - 1)  (ol 


(A.38) 


(O,  -1 


(A.39) 


(A.39): 


Equation  (A.37)  may  be  further  reduced  by  applying  Equations  (A.38)  and 


D’Azzo  and  Houpis,  p.  698  and  Laplace  Transform  15 
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m=- 


f  -^0)A  \ 

(  Q^dt  _  y 

f  -^laj  y 

CO, 

1 

+ 

K  t’t>d  J 

1  2  ; 

1  t^d  J 

1  2  J 

m 

m= 


0). 


(0„O), 


CO. 


n  y'^n'^d  J 


3l 

CO, 


{  A 


com. 


f  1 
e  ''  +e 


J  y^n^d  J 


dt\ 


1 


CO. 


f  \ 


\  CO  CO.,  I 

n  \  n  a  J 


<'  Q^d‘  _  Q-0>d<  \ 


■,/f^ 


f  ^ 


yf^nf^d  J 


f  Q(»dt  ^  ^ 


(A.40) 

(A.41) 

(A.42) 


Noting  the  identities  of  the  hyperbolie  sine  and  eosine  are^^; 
eosh(<»^t) : 


0),,t  ,  -COA 

e  ‘‘  +e  ‘‘ 


sinh(<»^t)  = 


(0)A  -0)A  \ 

e  ‘‘  -e  "  j 


(A.43) 

(A.44) 


Equation  (A.42)  may  be  simplified  by  appropriately  substituting  equations 
(A.43)  and  (A.44): 

r 

1 


fit)- 


co: 


(  qAo>J  \ 


yO^nO^dj 


^sinh(ry/)  +  (7?^)  oosh(<»^t) 


(A.45) 


Reealling  from  equation  (A.33): 

^co{t)'\  f  k.'^ 


y  Kd  y 


yJLj 


ir‘F(5) 


The  general  solution  for  angular  velocity  can  be  expressed  by  substituting 
equation  (A.45)  for  the  Inverse  Laplace  Transform  in  equation  (A.33): 


co{t)  = 


(KA] 

1 

1  JL  J 

_ 1 

{o}„co,)\ 

^^sinh(ry/)  +  (7^"-l)  cosh(<»^t) 


(A.46) 


To  develop  analytical  expressions  for  ^  and  (Dd,  the  earlier  definitions  of 
equation  (A.46)  must  be  revisited.  From  equations  (A.25)  and  (A.27)  we  know  that: 


2^«„  = 


V 


BL  +  RJ 
LJ 


y 


BR  + 


T-'-f 


y 


LJ 


J 


Douglas  F.  Riddle,  Calculus  and  Analytic  Geometry.  Alternate  Ed.,  Belmont:  Wadsworth,  1984, 

p.  406 
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Solving  for  ^  results  in; 


(  BL  +  RJ'^ 

2 

t  2JL  j 

(  BR  +  k^kf 

I  JL  ) 


(  BL  +  RJ^ 

1 

r  JL  ] 

t  2JL  j 

[BR  +  k^kf  j 

Solving  for  (Od  results  in: 


f  BR  +  k^kj^) 

f  BL  +  RJ'^ 

2 

^  JL  ^ 

-1 

1  J 

1  2JL  J 

^  BR  +  k^kf  J 

(  BL  +  RJ'^ 

2 

(  BR  +  k^k^'] 

t  2JL  J 

1  JL  ] 

co:, 


({Blf  +  2BRJL  +  {Rjf )  -  4JL  (bR  +  ) 

4(:Zf 


(A.47) 

(A.48) 

(A.49) 

(A.50) 

(A.51) 


0), 


If  BL-RJ^ 

2 

f  1 

(1  2JL  j 

1  JL  J 

(A.52) 


d.  Over-Damped  Angular  Velocity  Response 

The  angular  velocity  response  of  the  centrifuge  and  its  prime  mover  are 
greatly  dependent  upon  the  magnitude  of  If  ^>1  then  (Dd,  as  shown  by  equation  (A.27) 
will  be  a  real  number.  The  system’s  response,  as  governed  by  equation  (A.46),  will  be 
asymptotic  in  character,  or  ‘over-damped’.  The  over-damped  result  is: 


ty^>i(0  =  - 


KA 

\ 

1 

1 

JL 

a/  CO  CO. 

n  n  a 

KsK  1 


JL  (o: 


1- 


^  - 
CO, 


^  sinh  {o)A)  +  cosh  ( ) 

^  sinh  {o)j)  +  cosh  ( ) 


(A.53) 


(A.54) 


Recall  from  equation  (A.25)  that: 


col  = 


^  BR  +  k^kf 
LJ 


7 
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Substituting  for  ©n’^  yields; 

_nA  JL  L 

JL  {^BR  +  k^kj^^  0)^ 


\  1-— A®"'  ^ sinh  ( <»/ )  +  - 1  cosh  ( <»/ )  (A.55) 


co^>i{t)  =  -, — — r  1-  — 
^  {BR  +  k,k^)l  0)^ 


1-— e  ^sinh(<»^t)  +  -^^^  -lcosh(<»^t)  (A.56) 

0),  L  J 


Recall  from  equation  (A.  15)  that: 

CO  = — - 

(RB  +  k,kf) 

Therefore,  the  final  form  of  the  over-damped  response  equation  is: 

l--^e“^®'’'  ^sinh(ft»^t)-i-^^^  -lcosh(ft»/) 

Q)^  L 

e.  Critically-Damped  Angular  Velocity  Response 
If  ^=1,  then: 

Thus,  the  system’s  angular  velocity  response  will  be; 


(A.57) 


co.^it) 


1 

1 

1 

.^n 

(0„C0,  L 

1 

1 

1 

- 1 

Ct/ 

CO 

n 

n 

(A.58) 


(A.59) 


(A.60) 


KsK 

.  1 1  tn 

1 

1 

1 

JL 

iiili 

(O^  (0 

n 

n 

^sinh| 


cojJi-1 


- 1  cosh  uyt  - 1 


(A.61) 


co.^it)- 


KsK 

1 

1 

1 

JL 

co^ 

" 

CO 

n 

(A.62) 
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Using  substitution,  the  limit  of  the  hyperbolic  sine  term  was  found  to  be 
indeterminate  (the  value  found  was  zero  over  zero).  L’Hopital’s  Rule^i’^2  ^^vas  applied  to 
evaluate  the  limit’s  result: 


lim 


=  lim 


d 

d^ 

sinh 

1 - 1 

1 

<N 

d 

d^ 

1  1 

K) 

1 

1  1 

^ sinhi -1 
lim - ,  - 


lim^ 


( cosh 


d 


lim 


^sinh(ry„t7?^) 


=  cosh(0)t  =  t 


(A.63) 


(A.64) 


(A.65) 


Substituting  the  evaluated  result  of  the  hyperbolic  sine  term  from  equation 
(A.65),  equation  (A.62)  reduces  to: 


Vk. 


JL 


1  e 


(0„  co„ 


t  +  - 


1 

(0„ 


Equation  (A.66)  may  be  further  reduced  as  follows: 

V.X  1 


co,X)  = 


JL  (o:. 


-■l-e-^"”^Kt  +  l]] 


Recalling  the  value  of  (Dn’^  from  equation  (A.25): 


ty.  i(0  = 


VX 


JL 


JL  [BR  +  kJif] 


(BR  +  k^k^y 


(A.66) 


(A.67) 


(A.68) 


(A.69) 


Recalling  the  value  of  (Dss  from  equation  (A.  15)  results  in  the  final  form  of 
the  Critically-damped  Angular  Velocity  Response  Equation: 


D’Azzo  and  Houpis,  p.  130 
Riddle,  pp.  667-674 
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(A.70) 


(0  =  + 1]] 

f  Under-Damped  Angular  Velocity  Response 
When  ^<1,  then  (Da  beeomes  imaginary: 

(Od  =  (On  Vf  -1  =  =  j(On  (A.71) 

Substituting  equation  (A.71)  into  the  general  form  for  angular  veloeity 

yields: 


(0,<i{t)  =  0),, 


co^e  ^  " 


M 


(O,At)  =  C0,, 


\-e 


-^0)J 


^sinh cosh 


(A.72) 


J 


sinh 


7 


+  eosh 


(A.73) 


V  V 

The  imaginary  eomponent  may  be  eliminated  from  equation  (A.73)  by 
substitution  using  the  following  identities'^; 


•  ,  /  •  \  \  sinhlyh)  . 

sinh \ju)  =  7 sin (m )  ^ - =  sm (w ) 

j 

cosh(7M)  =  eos(M) 


u  =  cot 


The  result  is: 


(0,<At)  =  0),, 


-^0)nt 


1- 


Taking  advantage  of  the  Angle-Sum  Relation  of  Sine^'^: 

sin(M  +  ^)  =  sin  u  eos  (j)  -I-  eos  u  sin  (j) 

And  the  Quotient  Relation  of  Tangent^^; 


(A.74) 

(A.75) 

(A.76) 

(A.77) 

(A.78) 


William  H.  Beyer,  CRC  Standard  Mathematical  Tables.  28*  Ed.,  Boca  Raton:  CRC  Press,  p.  172 
Beyer,  p.  139 
Beyer,  p.  138 
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(A.79) 


tan  (p  =  ■ 


Then  the  like  terms  in  sine  and  eosine  may  be  equated  as  follows: 

cos^  =  <^ 


(f)  =  tan  ' 


(A.80) 

(A.81) 


(A.82) 


And  the  sine  and  eosine  terms  within  equation  (A.77)  may  be  redueed  as 


follows: 


sin =  +  (^l-^^Jeos (A.83) 

Reealling  equations  A.83,  A.84,  and  A.85: 

sin 1^1-^^  jt  +  (^)  =  cos((^)sm|ft»/^l-^^  j  +  sin [(j)) cos{o)J^j\-^^  j  (A.84) 

The  result  being  a  pair  of  equations  that  further  reduees  equation  A.80  and 
is  the  final  result  for  the  Under-damped  Angular  Veloeity  Response: 


=  1-^' 


sin  aU\-^  U  +  ip 


(A.85) 


(p  =  tan  ' 


(A.86) 


g.  Over-Damped  Current  Response 

Knowing  the  prime  mover’s  instantaneous  eurrent  requirements  is 
important  to  understanding  the  overall  power  requirements.  The  eurrent  to  voltage 
transfer  funetion  ean  be  obtained  through  the  substitution  of  equation  (A.7)  into  equation 
(A.10): 


l[s)  0)[s) 

V{s)  lty(5)llF(5) 


(A.87) 
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I{s)  _f  Js  +  b'^ 
v(s}  y  K  j 


{hA 

{jL) 


BL  +  RJ 
JL 


5  + 


^  BR  +  k^k^'^ 


JL 


(A.88) 


JJ 


zw 

V{s) 


JL  / 


f 

+ 

V  V 
f 


BL  +  RJ 
JL 


5  + 


^  BR  +  k^k^^ 
JL 


(A.89) 


V{s) 


^1} 

\Lj 


(  A 
I  J) 


'BL  +  RJ 
V  JL 


5  + 


V  y 


(A.90) 


Noting  the  same  eharaeteristie  equation  as  found  in  equation  (A.23),  the  ^ 
and  (Dn  terms  found  in  equation  (A.25)  may  again  be  substituted; 


M 

V[s) 


\Lj 


^  ^  B 

s  + 


V 


J 


+  {2^(u„)s  +  ((ul) 


(A.91) 


Reealling  the  definition  for  (Dd  found  in  equation  (A.27): 


M 

V{s) 


L 


5  +  - 


B 

J 


{s  +  {^0)„-0),)){s  +  {^0)„+0),)) 


(A.92) 


Onee  again,  we  assume  a  step  input  for  the  system  voltage  as  found  in 
equation  (A.  11).  This  results  in  a  step-response  transfer  funetion  for  the  eurrent  drawn  by 
the  system’s  prime  mover: 


L(.) 


L 


5-1- 


(s  +  {^a)„-a),))(s  +  {^a)„+a),)) 


(A.93) 


Taking  the  Inverse  Laplaee  Transform  of  equation  (A.93)  yields: 
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V  Ks  J 


L 

V  V 


s  +  — 

.  J . 


;(s  +  {<^co„-co,))(s  +  {<^co„+co,)) 


(A.94) 


L 


5  +  - 


B 

7 


^{s  +  {^0)„-0),)){s  +  {^0)„+0),)) 


(A.95) 


Referring  to  the  Inverse  Laplace  Transform  Tables^^  yields  the  following 


expressions: 


m 


F{s)  = 


[a-a)e 


s  +  a 


(5  +  a)(5  +  h)(5  +  c) 


^  {a-b)e  ^  {a-c)e 


(A.96) 


(A.97) 


{b-a){c-a)  {c-b){a-b)  {a-c){b-c) 

The  subtle  difference  between  equation  (A.97)  and  equation  (A.35)  should 
be  noted.  The  effects  of  friction  and  moment  of  inertia  are  much  more  profound  when 
predicting  current  than  angular  velocity. 

Finding  the  simplest  solution  of  the  Laplace  Domain  Component  of  the 
general  current  response  solution  follows: 

7  /  D  A  ^ 


/(0  =  ^‘ 


B 

s  +  — 

.  J . 


:{s  +  {<^a)„-a),]){s  + {<^0)^  +  0),]) 


(A.98) 


/W  = 


J 


.(A.99) 


f{t)= 


J 


2(0, 


2(0, 


({^(0„f-(0^,) 


.(A.  100) 


D’Azzo  and  Houpis,  p.  698,  Laplace  Transform  16 


86 


/(0  = 


/  \ 

{  B/ 

/  \ 

\ 

/  \ 

r  B/ 

/  \ 

\ 

/j 

* 

/j 

* 

^-{^(0„J(0^)t 

V 

V 

1 

/ 

_  J _ 

1 

\ 

1 

/ 

2(0, 


2(0, 


((^(0„f-(0^j) 


.(A.lOl) 


Noting  that: 


y  -  oj]  ^^col  -  col  - 1) 


(A.  102) 


(A.  103) 


02n 

Equation  (A.lOl)  may  be  reduced  by  applying  the  substitutions  noted  in 
equations  (A.  102)  and  (A.  103): 


/w=- 


B 

7 


(B/ 


0)^ 


-1 


-[40)„+(0j)t 


2V^ 


2V^ 


(A.  104) 


/(0  = 


5/ 
— - 

V,  "  ) 

e^‘‘‘ 

( 1 

( B/  '' 

K  "  J 

-<0.0 
e  ‘‘ 

1  2 

2  (A.  105) 

/0)=- 


B 

f  e  7 

-^(oj 

e  ^  " 

1  ) 
-(OA 

r 

V 

2 

2 

) 

(A.106) 


.g 


/(O- 


.6 


V  7ry„ 


Jco. 


n  J 


B 


B 


V  7<y„ 


Jco. 


n  J 


(A.  107) 


(y„ 


.8 

Jco„ 


-e 


-^0)J 


fit) 


B 


^  cosh{coj)  +  -^2L^smh{coj) 


Jco. 


co„ 


(A.  108) 
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Substituting  the  reduced  Inverse  Laplace  Transform  back  into  equation 


(A.95): 


B 

Jo)„ 


-e 


^\t) 


B 


V  L(0 

ss  n 

Factoring  out  the  common  coefficient  yields: 


(A.  109) 


K.  JLcol 


\-e 


(  r  JO), 
cosh(<»^t)  +  ^ 


B 


sinh(<»^t) 


(A.  110) 


The  leading  coefficient  can  be  reduced  by  using  values  previously  defined 
in  equations  (A.  17)  and  (A.25): 

J}  \(  TT  \  D 

(A.111) 


B  _| 

(  B  ^ 

JL  ^ 

_  B  _  4 

JLcol 

(JL, 

BR  +  k^kf  j 

BR  +  Kk^ 

Further  reducing  equation  (A.l  10)  to  the  following: 


^(0 


L. 


V  V 

55"  ss 


l-e 


cosh(<»^t)  + 


^  ^  Jco^ 
^  B 


sinh(<»^t) 


(A.  11 2) 


Therefore,  the  general  and  over-damped  solution  for  current  is: 


'■(0  =  4 


l-e 


cosh(<»^t)-i- 


_ _ ^ 


sinl 


h{coj) 


(A.  11 3) 


L  V  J] 

In  comparison,  the  general  form  of  the  current  response  (equation  (A.l  13)) 
contrasts  with  the  general  form  of  the  angular  velocity  response  (equation  (A.46))  only  in 
the  coefficient  of  the  hyperbolic  sine  term. 

Equation  (A.46)  is  repeated  here  for  the  purpose  of  comparison: 


88 


fV  k  ^ 

«(0=  ^ 

\  JL  J 


(0„ 


\0)n0)d  J 


^sinh(6;/)  +  (V^"-l)  cosh(<»^^) 


h.  Critically-Damped  Current  Response 

Noting  the  similarities  between  the  general  form  of  the  eurrent  response 
(equation  (A.l  13))  and  the  general  form  of  the  angular  velocity  response  (equation  (A.46) 
),  the  critically-damped  current  response  can  be  deduced  from  the  critically-damped 
angular  velocity  response.  Recalling  from  equation  (A. 70): 


Therefore: 


(0  =  4 


( 

\-e  " 

Ti 

\ 

cot  +  \ 

\ 

B  j 

n 

J_ 

(A.  11 4) 


L  Under-Damped  Current  Response 

As  with  the  general,  over-damped,  and  critically  damped  cases;  the  under¬ 
damped  current  response  can  be  deduced  from  the  under  damped  angular  velocity 
response.  Recalling  from  equation  (A. 77): 


1- 


Realigning  the  coefficients  in  order  to  provide  easier  deduction  yields: 


\-e 


-^0}J 


cos 


^  '-I-  ,  ^  sin 


(A.  11 5) 


Note  the  coefficient  of  the  sine  term  in  equation  (A.  11 5)  is  the  inverse  of 
the  phase  angle  terms  in  equation  (A. 86). 


(f)  =  tan  ' 


I 


Through  deduction,  the  under-damped  current  response  is: 


89 


(A.  11 6) 


k<y  (0  =  hs 


\-e 


-^(Ont 


COS 


+ 


r  B 


sm 


Expressed  in  terms  of  phase  angle; 


-^0)J 


rSm 


(/>  =  tan 


1  ^/W 


r  B 


j.  Collected  Results 

Angular  Veloeity  Transfer  Funetion,  equation  (A.22): 


co{s) 


K. 


BR  +  k^kj^  j 


\BR  +  k,kf) 

LJ 

f  BL  +  RJ'^ 

5  + 

(  BR  +  k^k^'\ 

V  LJ  J 

1  LJ  JJ 

(-] 

UJ 


Steady  State  Angular  Veloeity  Response,  equation  (A.  15): 


0),.  = 


(RB  +  k,kf) 


K. 


Natural  Frequeney,  equation  (A.25): 


^■.2 

= 


f  BR  +  k^k^'^ 
LJ 


Damped  Natural  Frequeney,  equation  (A.27): 

Damping  Ratio,  equation  (A.48): 


(  BL  +  RJ^ 

1 

JL  ■] 

t  2JL  J 

BR  +  k^kf  ^ 

(A.117) 

(A.  11 8) 
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Over-Damped  Angular  Velocity  Response,  equation  (A. 57): 


CO, 


^ sinh {(Oj)  +  -1  cosh ( ) 


Critically-Damped  Angular  Velocity  Response,  equation  (A.70): 


0) 


Under-Damped  Angular  Velocity  Response,  equations  (A.85)  and  (A.86): 


J 


0  =  tan  ' 


Jhl 

I 


Current  Transfer  Function,  equation  (A. 90): 


M 

V{s) 


(  1  ^ 


B 

s  +  — 

\  J  J 


5^-1- 


BL  +  RJ 

jzT 


^  f  BR  +  k^kf'^ 


5-1- 


JL 


Current  Steady  State  Response,  equation  (A.l  1 1): 

i  =V  - - - 

”  BR  +  k^kf 

Over-Damped  Current  Response,  equation  (A.l  13): 

r  f  J...  \  \ 


C>i 


(0^ 


\-e 


cosh(<»^t)-l- 


L! _ B_j_ 


smh(<»^t) 


Critically-Damped  Current  Response,  equation  (A.l  14): 


/ 

1  -(k>J 

\-e  " 

Ti 

\ 

(Ot  +  \ 

\ 

B  J 

n 

)  _ 

91 


Under-Damped  Current  Response,  equation  (A.l  17)  and  (A.l  18): 


C<i(0  =  4 


-^0)J 


rSin 


(j)  =  tan 


,  ^/w 


''t  ■/<». 

^  B 


3,  State-Space  Formulation  of  the  Centrifuge  External  Dynamics 
Problem 

A  second  means  of  solving  the  relationship  for  current  and  angular  velocity  was 
needed  in  order  to  ascertain  the  validity  of  the  Analytical  Model.  Recalling  from 
equations  (A.l)  and  (A.2): 


Ri  +  L 


di 

dt 


=  v-kj- 


dt 


kj  =  J 


d^e  „  de 

— r  +  B  — 

dt  dt 


resolution: 


Rewriting  in  terms  of  to: 


Ri  +  L  —  =  V-k,0) 
dt  ^ 


rdO)  n  ,  . 

J - h  Bo)  =  kj 

dt 


Solving  for  the  derivatives  yields: 


di 

dt 


( R  .  ^ 

—  1  +  —0) 

U  L  ) 


V 

+  — 
L 


do) 

dt 


J  J 


(A.  11 9) 
(A.  120) 


(A.121) 
(A.  122) 


Rewriting  equations  (A.121)  and  (A.  122)  in  matrix  form  for  simultaneous 


d  i 

i 

(R/L) 

i 

4- 

(l/i) 

dt  c 

0 

(KIJ) 

-(BIJ)  . 

0) 

0 

(A.  123) 
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Equation  (A.  123)  may  be  represented  by  the  state-spaee  formulation: 


X= AX+BU 


Where: 


-  ,  I 

X 

'a\ 

^ R')  f 


V  ^  J 


\Lj 


\J  J 

B 


V  Jj 


0 


And: 


U  =  E 


At 


Substituting  equation  (A.  129)  into  equation  (A.  124)  yields: 

l  +  BU 


At 


Simplifying  equation  (A.  130): 


X„„-X,=AA( 


^X„^+X,^ 


J 


+  BUAt 


X,,..=AAt 


X,,..-AAt 


f  V  V  ^  _ 

‘  l  +  BUAt  +  X, 


V 


J 


fx  .  ^ 

fx  ^ 

=  AAt 

1  2  ; 

l2  j 

AAt 

I  ^ 


X, 


AAt 

1  + - 

V 


+  BUAt  +  X, 


X,+BUAt 


J 


X.,..  = 


(  AAt^ 

-1  ■ 

Y  AAt^ 

I 

1  + 

X,+BUAt 

1  2  ; 

.1  2  ) 

t 

(A.  124) 

(A.  125) 

(A.  126) 

(A.  127) 
(A.  128) 

(A.  129) 

(A.  130) 

(A.131) 

(A.  132) 

(A.133) 
(A.  134) 

(A.135) 
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The  two  terms  in  equation  (A.  135)  containing  the  identity  matrix  must  be 
resolved  before  the  determinant  can  be  found.  First: 


(A.  136) 


(A.  137) 


(A.  13  8) 


(A.  139) 


(A.  140) 


(A.141) 


The  inverse  of  a  matrix  can  be  found  using  the  following  rule^^: 


M 


2x2 


a  b 
c  d 


Dennis  G.  Zill,  A  First  Course  in  Differential  Equations  with  Applications.  4*  Ed.,  Boston:  PWS- 
Kent,  1989,  p.  384 
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M 


-1 

2x2 


d  -b 
-c  a 
[ad) -{be) 


Therefore,  the  seeond  term  is: 


AAt 


(,  BAt^ 

1  + - 

^  2J 


^  k^At^ 


I  2J  ) 


V 

( 

1- 

l  2Ty 


J 

RAt'^ 


I- 


"l_  — 
V  2L 


BAt 

1  + 


^  f kfAt^k 


AAt 


f 

V 


2J 


1  + 


BAt 


2J 

(  kAt\ 


V  yv 


^  kfAt^ 


k^At 


2J 


2L  j 
RAt'^ 


2L 


.  RAt  BAt  BRAt 
1 - +  - 


2L  2J 


4JL 


kj^k^At 

4JL 


2  A 


(A.  142) 


(A.  143) 


Recalling  the  terms  defined  in  equations  (A.  125),  (A.  126),  (A.  127),  (A.  128), 
(A.  138),  and  (A.  143)  yields  the  following  when  substituted  into  equation  (A.  135): 


X, 


1  + - 

t  2J 

'KAt'^ 

I  2J  j 


''  k^At  ^ 


2L  J 
RAt'^ 


2L 
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I  2T  j 
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I  2J  j 


''  kfAt^ 

2L  J 
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1- 


2J 


+ 


UJ 

0 


VAt 


,  RAt  BAt  BRAt 
1 - + 


X.,„  = 


BAt 

1  + - 

2J 

f  kAt\ 


2J 


2L  2J 


^  kfAA 


1- 


2L  J 

raA 


2L 


i  + 


k^kA^ 

4JL  4JL 
iRAt  kfCOAt 


2  A 


2L 
ikAt 


2J 


+  0)- 


2L 

BcoAt 


2J 


+ 


'vaA 

I  L  j 


t+At 


^  RAt  BAt  BRAt 
2L  2J  4JL 


k^kA^ 

4JL 


2  y 


(A.  144) 


(A.  145) 
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<1 

i,RAt  kfCoAt  VAt 
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V 

2L  J 
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V 
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(  k,At 


2L  2J  4JL 
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(A.  146) 


4JL 

(_kf^](  AAt 

2J 


2J 


BcoAA 
2J  j 

A 


2J  ) 


^+A? 


,  RAt  BAt  BRAt^ 
1 - + 


2L  2J 


4JL 


k^At 

4JL 


2  A 


(A.  147) 
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1  + - 

^  2J 
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\RAt  kfCoAt  VAt 
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2L 


2L 
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L 


ifk^A^  kfCoA^  BkfCoAt 


2  A 
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4JL 


2L 


4JL 


J 


2  T  xr  »  ,2  ^  / 


LkA^  i.kA^t^  kfk^wAt  ky  At 

^  2J  ^  4JL  JjL  ^  2JL  J 


+ 


,  7?A^V  i,k,At  BcoAA 
1 - —  +  0).  -  ‘ 


V 


2L 


JK 


2J 


2J 


J 


^ ^  RAt  ^  BAt  BRAt^  k^A^^^ 
^  2L  2J  4JL  4JL  j 


(A.  148) 


The  simultaneous  solution  of  and  0)^^^  yields  the  final  form  of  eaeh 
equation. 


1  + 


BAA( .  i,RAt  k.coAt  VAA  ( i.kA^At^  kfCO.At  Bk,co,At^^ 


2J 


I*  “I"  ■ 


^t+At 


2L 


V  j'KAy^L. 


-+- 


L 


4JL 


-  + 


Ax 


V  j'KAy^L. 

4JL 


^ ^  RAt  ^  BAt  BRAt^  k^A^^^ 
2L  2J  4JL  4JL  j 


(A.  149) 


i,kAt  _l_  i,kA^^^  kykr^A^^  _l_  k^V At^ 


A+Ai 


_  V 


2J 


4JL 


4JL 


2JL 


+ 


1- 


J 


V 


RAt 

AI 


\r 


j\ 


i,k,At  BcoAt 

— —  +  0), - — 

2J 


\ 

2J  y 


RAt  BAt  BRAt^  kXA^^^ 

1 - + - kX - 

2L  2J  4JL  4JL  j 


(A.  150) 
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The  numerical  results  of  the  analytical  model  (Equations  (A.  15),  (A.  17),  (A.25), 
(A.27),  (A.48),  (A.57),  (A.70),  (A.85),  (A.90),  (A.lll),  (A.113),  (A.114),  (A.117),  and 
(A.118))  and  the  state  space  model  (Equations  (A.15),  (A.17),  and  (A.123))  have  been 
compared  and  are  equivalent. 
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B,  APPENDIX  B 

1,  Modeling  the  Internal  Dynamics  of  a  Centrifuge 

This  section  discusses  the  mathematical  rigor  behind  the  dynamics  of  a  rotational 


environment. 

Equation  Section  (Next) 


Analytical  Term 

Meaning 

Units 

iR 

Centrifuge  Radial  Axis 

Unitless 

iv 

Centrifuge  Angular  Axis 

Unitless 

ik 

Centrifuge  Axial  Axis 

Unitless 

R 

Linear  Position 

m 

V 

Linear  Velocity 

m/s 

e 

Angular  Position 

rad 

0) 

Angular  Velocity 

rad/s 

g. 

Average  Earth  Gravity  at  Sea  Level,  9.81  m/s^ 

m/s^ 

gc 

Centrifugal  Acceleration 

m/s^ 

G 

“G”  Level 

Unitless 

F 

Force 

N 

Table  15:  Centrifgue  Internal  Terms 


2,  Velocity 

The  velocity  of  a  subject  (as  shown  in  Figure  61)  within  a  centrifuge  may  be 
expressed  in  terms  relative  to  the  external  world  as  follows: 

r.,=^+(«xR)  (B.i) 

3,  Acceleration 

Similarly,  the  acceleration  of  a  subject  within  a  centrifuge  may  be  expressed  in 
terms  relative  to  the  external  world  as  follows: 

/  dW ext  \  /T-> 

— - = - +  tax  Vert  (B.2) 

m  dt  ^ 
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Substituting  the  equation  (B.l)  into  equation  (B.2)  yields: 


y  F 

/  j  ext  _ 


d 


^dK  /- 


dt 


-  + 


(ojxr) 


m 


dt 


-  + 


JR  /—  — \ 

ox 

- -1-  ojxR 

V 

dt  ^  ’ 

V  Jj 

(B.3) 


Expanding  equation  (B.3)  and  collecting  terms  results  in  the  following: 


y  F 

/  V  ext  _ 


{  f  J _ 


m 


JR 


-  + 


Je>  — 

- xR 

dt 


W 


+ 


JJ 


-  JR 

2e)X - 

dt 


+ 


(ox^wxr)!  (B.4) 


Figure  61:  Centrifuge  Internal  Terms 


a.  Special  Condition:  Micro-Gravity,  No  Thrust 
In  a  micro-gravity  environment,  there  will  normally  be  no  external  forces 
acting  on  the  centrifuge.  Therefore: 


m 


(B.5) 
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b.  Special  Condition:  Centrifuge  Steady  State 

Once  the  eentrifuge  is  “up  to  speed”  angular  veloeity  will  remain  eonstant. 

Therefore: 


(B.6) 


c.  Application  of  Special  Conditions 

Under  the  speeial  eonditions  detailed  in  equations  (B.5)  and  (B.6), 
equation  (B.4)  may  be  redueed  as  follows: 


I  (  //  R  I 

+  (0xR)  +  2e>x -  +(e>x(e>xR)j 

dt  , 


(B.7) 


dt^ 


+ 


(B.8) 


d^K 

df 


JR^ 

X - 

dt  ^ 


(B.9) 


4,  Apparent  Forces 

Equation  (B.9)  shows  three  groups  of  terms.  The  first  of  which  is  Apparent 
Aeeeleration: 

J'R 

The  seeond  term  is  Centrifugal  Foree: 

ox^Rxc)) 

The  third  term  is  Coriolis  Foree: 


a.  Centrifugal  Force 

The  axial  eomponents  of  centrifugal  foree  may  be  developed  by  utilizing 
Vector  Analysis: 
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Rxe> : 


Ir  Iv  Ik 

i?  0  0 

0  0  <» 


-Rcoi, 


(B.IO) 


OJX 


(Rxoj)  : 


Ir  Iv  Ik 

Q  Q  CO 

0  -Rco  0 


Ro/iR 


(B.ll) 


The  results  of  equation  (B.ll)  may  be  validated  through  eomparison  with  two 
referenees^^’^^  and  by  reealling  equation  (B.5): 


centrifugal 


m 


(B.12) 


b.  Coriolis  Force 

As  with  centrifugal  force,  the  axial  components  of  Coriolis  Force  may  be 
developed  by  utilizing  Vector  Analysis: 


2e>x- 


dt 


=  -2 


Ir 

0 

C  r7D  \ 


dt 


Iv 

0 

C  r7D  \ 


dRj,]  dR^]  ( dR, 


dt 


Ik 
CO 

C  riD  \ 


_ k 

dt  j 


=  2 


f 

^  dR,^ 

f  dR^^ 

_  A 

CO 

iR  -  CO 

iv 

V 

^  dt  J 

dt  j 

y 

(B.13) 


The  results  of  equation  (B.13)  may  be  validated  through  the  use  of  two 
references®0’6i  and  by  recalling  equations  (B.5)  and  (B.6): 


F 


f 


Coriolis 

m 


dR- 


CO 

V  dt 


V  ■  dR;,- 

iR  -CO — —I 


dt 


(B.14) 


Note  that  equation  (B.14)  shows  that  there  are  two  components  of  Coriolis 
Force  in  a  centrifuge.  The  first  is  the  result  of  a  linear  change  in  angular  position  (for 
example,  walking  a  distance  across  the  floor  into  or  away  from  the  rotation  of  the 


J.L.  Meriam  and  L.G.  Kraige,  Engineering  Mechanics  Vof  2:  Dynamics.  2”“*  Edition,  New  York: 
John  Wiley  and  Sons,  1986,  pp.  223-224  and  350 

Dare  A.  Wells  and  Harold  S.  Slusher,  Schaum’s  Outline  of  Theory  and  Problems  of  Physics  for 
Engineering  and  Science.  New  York:  McGraw-Hill,  1983,  p.  39 

Meriam  and  Kraige,  Dynamics,  p.  350 

61  Wells  and  Slusher,  p.  64 
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centrifuge).  This  foree  appears  as  either  an  inerease  or  deerease  in  apparent  gravity,  as 
follows: 

F  ///? 

Coriolis  .horizontal  r\  "T  /y-j 

- =  2o) — -iR  (B.15) 

m  dt 

The  second  component  of  Coriolis  Foree  is  the  result  of  a  linear  change  in 
radial  position  (examples  inelude:  sitting,  standing,  or  elimbing  a  ladder).  This  force 
appears  as  a  shear  foree  that  oecurs  with  or  against  the  rotation  of  the  centrifuge,  as 
follows: 

F  ///? 

Coriolis , vertical  Cil\.n  “ 

- =  -2(0 — -iv  (B.16) 

m  dt 


5,  Vertical  Coriolis  Versus  Centrifgue  Radius 

Focusing  on  the  disorienting  effects  of  vertical  Coriolis  on  human  equilibrium,  a 
relationship  between  vertical  Coriolis  and  centrifuge  radius  should  be  developed. 
Redefining  centrifugal  acceleration  in  terms  of  G-Level  yields: 

m 

_  ^centrifgual  _  R(0  ^ 

o  centrifugal  v  ’  ) 

go  go 

Therefore,  to  find  the  proportionality  of  vertical  Coriolis  to  centrifugal  radius,  the 
following  expression  must  be  developed: 


Coriolis  ..vertical 


S  centrifugal 


- T^r^lv  = - ^ - -Iv 


(B.19) 


Redefining  centrifugal  acceleration  in  terms  of  angular  velocity  yields: 


Thus,  equation  (B.19)  may  be  redefined  as: 


(B.20) 
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f  F  ^ 

Coriolis , vertical 


m 


j  _ 


2. 

\  cit  J 


S  centrifugal 


R. 


S  centrifugal  S  0 

R 


F. 


Coriolis  .vertical 


m 


J  _ 


=  2 


S  centrifugal 

^  F 


( 

fgo' 

^  R 

y  dt  J 

Vrj 

S  centrifugal  S  0 

V 


Coriolis , vertical 

m 


J  _ 


=  2 


S  centrifugal 


dR^ 

dt 


go 


centrifugal 


(B.21) 


(B.22) 


(B.23) 


Therefore,  for  a  given  “G-Level,”  the  vertieal  Coriolis  is  inversely  proportional  to 
the  square  root  of  the  centrifuge  radius. 


6,  Gravity  Gradient 

Developing  a  relationship  that  defines  gravity  gradient  is  useful  in  scaling 
experiments  so  that  their  applicability  to  a  human’s  experience  of  a  centrifugal 
environment  is  best  reproduced.  Recalling  from  equation  (B.18): 

_  Ro)^ 

g  centrifugal 

go 


The  radial  change  in  apparent  gravity,  gravity  gradient  is  the  derivative  of 
equation  (B.18)  with  relation  to  radial  position: 


dR^^centrifugal) 

(B.24) 

centrifugal  ^  ^ 

dR  go  go 

(B.25) 

To  find  the  proportion  of  gravity  gradient  to  centrifugal  radius,  the  following 
expression  must  be  developed: 


S centrifugal  _ 

S  centrifugal 


Uq  J 

Ro/^  R 


\  go  ) 


(B.26) 
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Therefore,  for  a  given  “G-Level,”  the  gravity  gradient  is  inversely  proportional  to 
the  centrifuge  radius. 
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c. 


APPENDIX  C 


1,  Modeling  the  Human  Ability  to  Sense  Rotation 

This  section  discusses  the  mathematical  rigor  behind  the  human  disorientation 


model. 

Equation  Section  (Next) 


Analytical  Term 

Meaning 

Units 

b 

Torque  per  unit  of  relative  angular  velocity62 

1/s 

F, 

Force  on  the  head 

N 

iR.v^k 

Centrifuge  Radial,  Angular,  and  Axial  Axes 

Unitless 

A 

Mass  Moment  of  Inertia  of  the  head 

kg-m^ 

k 

Torque  per  unit  of  relative  angular  displacement 

1/s 

R 

Finear  Position 

M 

V 

Finear  Velocity 

m/s 

e 

Angular  Position 

rad 

Torque  on  the  head 

N-m 

Cupula  Angular  Position 

rad 

Cupula  Angular  Velocity 

rad/s 

V 

Head  Angular  Acceleration 

rad 

CO 

Centrifuge  Angular  Velocity 

rad/s 

Table  16:  Human  Perception  Terms 


2,  Coriolis  in  a  Centrifuge 

The  somatogravic  illusion  (shown  in  Figure  62)  of  level  ground  generated  by  the 
centrifuge  is  disrupted  whenever  an  Astronaut  moves  as  shown  in  Appendix  B,  equation 
(B.9): 


/—  /—  — 

f  -  t/R^ 

(OX  Rxo)  - 

2(ox - 

V  V  // 

dt  ) 

Milsum,  p.  187 
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Figure  62:  Human  Perception  of  the  Centrifugal  Environment 

As  the  Astronaut  Moves,  a  new  Somatogravie  Illusion  is  generated  due  to  the 
application  of  new  force  vectors.  The  added  force  vectors  result  from  the  Coriolis  Term 
in  equation  (B.9)  and  are  detailed  in  equation  (B.14): 


f 


Coriolis 

m 


dR~ 


(0 

V  dt 


V  ■  dR,,- 

iR  -CO — —I 


dt 


Solving  equation  (B.14)  in  terms  of  Force  and  expressing  position  change  in 
terms  of  velocity  yields: 


Coriolis 


■■  ImCOVJr  -  ImcoVJv 


(C.l) 


3,  Describing  the  Semi-Circular  Canals  in  Terms  of  Head  Rotation 

The  relationship  between  the  orientation  of  the  head  in  space  and  the  semi¬ 
circular  canals  is  given  by  the  following  second  order  equation  based  on  Newton’s 
Second  Law®^; 

b[y/-d)  +  k[y/-0)  =  je  (C.2) 


Milsum,  p.  187 
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Cupula  Deflection  ((|))  per  unit  relative  angular  displacement  between  the  skull 
and  endolymph  fluid  is  defined  as64; 

(l)  =  a{y/-0)  (C.3) 

The  coefficient,  a,  is  the  area  ratio  between  the  canal  and  ampulla,  however  this 
model  is  based  on  a  general  case  and  thus: 

a  =  1  (C.4) 

(l)  =  y/-0  (C.5) 

Rewriting  equation  (C.2)  as  a  consequence  of  equation  (C.5)  yields: 

J(p  +  b^  +  k(j)  =  Jifr  (C.6) 


4,  Torque  Imposed  on  the  Endolymph  hy  Coriolis 

The  endolymph  in  the  semi-circular  canal  resides  in  what  is  essentially  a  thin 
tube.  It’s  Mass  Moment  of  Inertia  (Jh)^^  is  best  described  as: 

J,=m,Rl  (C.7) 

Substituting  equation  (C.7)  into  the  right  side  of  equation  (C.6)  yields: 

Jipr  =  m^R^i^  =  Tf,  (C.8) 

Noting  the  definition  of  torque  (x)®^: 

t  =  FxR  (C.9) 

Determining  the  value  of  x  requires  that  the  vector  product  of  Coriolis  and 
distance  from  the  center  of  rotation  be  found: 


Ir  Iv 

2m^coV^  -2mf,coVj^ 
R  0 


:  Intf^coVj^Rik 


(C.IO) 


It  is  worth  noting  that  Xh  due  to  Coriolis  is  only  felt  in  the  axial  plane. 

Th  ik  = 


(C.ll) 


64Milsum,  p.  187 

Meriam  and  Kraige,  Dynamics,  p.  630 

Meriam  and  Kraige,  Engineering  Mechanics  Vof  1 :  Statics.  2"^  Edition,  New  York:  John  Wiley  and 
Sons,  1986,  p.  30 
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5.  Quantitative  Analysis  of  Angular  Acceleration  Inputs  Resulting  from 
Coriolis 

To  determine,  through  quantitative  analysis,  the  angular  acceleration  input  from 
the  semi-circular  canals  as  a  result  of  Coriolis,  the  solution  to  (j)  and  Cupula  Veloctiy 
must  be  found.  Recalling  equation  (C.8)  and  substituting  the  determinant  of  equation 
(C.IO)  yields  the  following: 


='^k-  '^■m^coV^R 

(C.12) 

Solving  equation  (C.12) 

in  terms  of  ifr : 

2m,0)V„R  2a>V„ 

y/  =  — - — ^  = - - 

m,R^  R 

(C.13) 

Substituting  the  result  of  equation  (C.13)  into  the  right  side  of  equation  (C.6): 

t2  ,, 

J(j)  +  b(l)  +  k(l)  =  J - — 

R 

(C.14) 

2  b  :  k  2cdV„ 

(/>+  ^ 

J  J  K 

(C.15) 

Allowing: 

II 

(C.16) 

By  substitution: 

b  k  2o)V„ 

X,  -I-— X,  -I-— X  = - - 

J  J  R 

(C.17) 

Noting  that: 

Xj  =X2 

(C.18) 

2o)Vn  b  k 

X'y  =  - - - X'y - X, 

R  J  J 

(C.19) 

By  letting: 


no 


X 

0 

V  Jj 


X, 


1 

V  Jj 


U: 


0 

2o)V, 


\\  R  j 

Equation  (C.20)  can  be  simplified  as  follows: 

X= AX+U 

Noting  that: 

,  -x 


At 


Equation  (C.24)  ean  be  expanded  as  follows: 


X,^,  -X 


At 


A 


X,^„  +  X 


+  U 


Solving  for  X^^^,  yields: 


X,,,,-X,=AAt 


'  +  U 


X...,=AAt 


V  V 

^  l  +  u  +  x. 


Removing  the  X^^^,  from  the  right  side  of  equation  (C.28): 


X.,.,  = 


AAtX,,,,  +  AAtX, 


^X.  ^ 


X.,..-AAt 


t+M 


I- 


AAt 


V  2  y 

X,^„ 


fx  ^ 

=  AAt  ^  +U  +  X, 

U  j 


1  + 


AAt 


X,  +  UAt 


X.,..  = 


^  AAt 

”V 

,  AAt 

I 

2  ) 

\ 

1  + 

2 

X^  +  UAt 

y 
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(C.21) 

(C.22) 

(C.23) 

(C.24) 

(C.25) 

(C.26) 

(C.27) 

(C.28) 

(C.29) 

(C.30) 

(C.31) 


(C.32) 


The  two  terms  in  equation  (C.32)  containing  the  identity  matrix  must  be  resolved 
before  the  determinant  can  be  found.  First: 


1 

0 

At 

0 

1 

I-l - 

2  J 

0 

1 

H - 

2 

1 

-m 

^  AAt 

T  1 

1  0 

1 

0 

1 

2 

'  0  1 

“T 

-(*%)  -('’%) 

_  AAtl 

1  + -  = 

1 

2  J 

Second,  the  inverted  term: 


^  AAt 

1 

0 

At 

0 

1 

I - 

2  J 

0 

1 

T 

1 

-{%) 

j  AAt 

1  0 

0  A^ 

2 

'  0  1 

-("%)  -(’’Vij) 

The  inverse  of  a  matrix  can  be  found  using  the  following  rule®^: 

\a  b\ 


M 


2x2 


d 


(C.33) 

(C.34) 

(C.35) 

(C.36) 

(C.37) 

(C.38) 


d  -b 


[ad) -{be) 
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Therefore,  the  seeond  term  is: 


Mt  i  [At 


bAt  ^  f  At  V  kAt 


2J  J  I  2J  Jl  2J 


(C.39) 


bAt [At 


.  bAt  kAt^ 

1  + - + - ^ 

2J 


(C.40) 


Recalling  the  terms  defined  in  equations  (C.16),  (C.21),  (C.23),  (C.32),  (C.35), 
and  (C.40)  yields  the  following: 

Ir,  bAA  ( At3| 


l+f  [-(*%)  [— ]l 


k(l>At/ 

njy^  nj 


.  bAt  kAy 

1  + + ^ 

2J  AA 


+  I  2coVg^At 


(C.42) 


k(l)At 


.  bAt  kAt^ 

1  + - + - ^ 

2J  AA 


b^A^ 


2coVAt 


(C.43) 
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bAt^ 

1  + - 

V  2J 


<!>,+ 


(l),At 


+ 


k(/),At  ^At  b^^At  2o)Vj^At 
4J  2  4J  2R 


<1 

1 

(  ,  ^,A^^ 

+  - — 

-I- 

1 

It 

y  2  ; 

V 

k^,At  ^  ^  b^^At  2coVj^At 


2J 


2J 


R 


t+M 


.  bAt  kAt^ 

1  + - + - ^ 

2J  4J^ 


(C.44) 


The  simultaneous  solution  of  and  yields  the  final  form  of  eaeh 
equation. 


= 


l-i- 

f  ,  ^,At^ 

-I- 

1  2J) 

r  2 ; 

V 

k^,At  ^At  b^^At  2o)Vi^At 

4J  2  4J  2R 


2  A 
) 


,  bAt  kAt^ 
1  + - + 


(C.45) 


(  kAA 

(  ,  ^,At^ 

-1- 

1  2jJ 

1  2  ; 

V 

2J  4J 
k^^At 


2J 


■  b^.At  2(oV„At 
+  ^ 


2J 


R 


.  bAt  kAt^ 
1  + - +  - 


(C.46) 


2j  4r 


6,  Analytical  Check  Case 

An  analytical  check  case  was  performed  using  Mathematica  version  4.^^  The 
numerical  results  of  the  Mathematica  model  and  the  state-space  model  (Equations  (C.45) 
and  (C.46))  were  compared  and  found  to  be  equivalent. 
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D, 


APPENDIX  D 


1,  Supporting  Subroutines 

Equation  Section  (Next)Several  LabView  Programs  were  written  to  simplify  the 
code  used  in  the  main  computer  models.  The  details  behind  each  supporting  sub-routine 
are  listed  in  the  following  sections. 

2,  Unit  Conversion  Subroutines 

The  natural  unit  for  angular  position  is  radians.  It  is  difficult  to  visualize  the 
6.283  radians  that  make  up  a  complete  circle.  Degrees,  Hertz,  and  RPMs  however,  are 
much  easier  for  the  user  to  conceptualize  and  communicate  to  a  target  audience. 
Therefore,  an  automated  means  of  performing  the  repetitious  conversion  process  from  a 
human  measure  (degrees)  to  a  natural  measure  (radians)  is  needed. 


Analytical  Term 

Meaning 

Units 

b 

Torque  per  unit  of  relative  angular  velocity®^ 

1/s 

Fu 

Force  on  the  head 

N 

iR.v.k 

Centrifuge  Radial,  Angular,  and  Axial  Axes 

Unitless 

Mass  Moment  of  Inertia  of  the  head 

kg-m^ 

k 

Torque  per  unit  of  relative  angular  displacement 

1/s 

R 

Linear  Position 

M 

V 

Linear  Velocity 

m/s 

6 

Angular  Position 

rad 

h 

Torque  on  the  head 

N-m 

Cupula  Angular  Position 

rad 

Cupula  Angular  Velocity 

rad/s 

¥ 

Head  Angular  Acceleration 

rad 

CO 

Centrifuge  Angular  Velocity 

rad/s 

Table  17:  Terms  Used  in  tbe  Unit  Conversion  Subroutines 


Milsum,  p.  187 
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a. 


Degrees  to  Radians  Converter 


This  VI  Converts  An  Angular  Position  Input 

From  Degrees  To  Radians  And  Vice-Versa 

Angle  Input  deg 

Angle  Output  rad 

0.000 

0.000 

Angle  Input  rad 

Angle  Output  deg 

0.000 

0.000 

Figure  63:  Degrees  to  Radians  Converter  Front  Panel 


The  front  panel  of  the  Degrees  to  Radians  Converter,  shown  in  Figure  24, 
is  not  designed  for  direet  user  manipulation.  Rather,  this  VI  is  designed  to  be  a 
subordinate  VI  invoked  whenever  a  routine  eonversion  between  degrees  and  radians  is 
required.  The  Angle  Input  Controls,  are  linked  from  an  outside  souree.  The  necessary 
conversion  computations  are  performed  with  the  results  passed  to  the  Angle  Output 
Indicators. 


0 


Figure  64:  Line  Diagram  of  the  Process  to  Convert  Degrees  to  Radians 
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Equation  (D.l)  was  used  to  perform  the  ealculation  shown  in  Figure  64 

0 

0  =27t^^  (D.l) 

360 


Equation  (D.l)  was  solved  for  degrees  resulting  in  equation  whieh  was 
used  to  perform  the  calculation  shown  in  Figure  65. 


360 

iTt 


0 


(D.2) 


b.  Frequency  to  Omega  2-Way  Converter 

The  front  panel  of  the  Frequency  to  Omega  Converter,  shown  in  Figure 
66,  is  not  designed  for  direct  user  manipulation.  Rather,  this  VI  is  designed  to  be  a 
subordinate  VI  invoked  whenever  a  routine  conversion  between  hertz  and  radians  per 
second  is  required.  The  Angular  Velocity  Input  Controls,  are  linked  from  an  outside 
source.  The  necessary  conversion  computations  are  performed  with  the  results  passed  to 
the  Angular  Velocity  Output  Indicators. 

Beyer,  p.  152 
Beyer,  p.  152 
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This  VI  Converts  An  Angular  Velocity  Input  From 

Hertz  To  Radians  Per  Second  And  Vice-Versa 

Angular  Velocity  Input  hz 

Angular  Velocity  Output  rad/sec 

^0.00 

0.00 

Angular  Velocity  Input  rad/sec 

Angular  Velocity  Output  hz 

^0.00 

0.00 

Figure  66:  Frequency  to  Omega  Converter  Front  Panel 


Figure  67:  Line  Diagram  of  Process  to  Convert  Hertz  to  Radians  Per  Second 


The  formula  used  in  Figure  67,  Equation  (D.3),  is  derived  from  the  fact 
that  there  are  2%  radians  in  a  complete  circle^^^  and  one  hertz  represents  one  revolution 
per  second. 


0)  =  2^f 


(D.3) 


Beyer,  p.  152 
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Solving  equation  (D.3)  for  /  yields  equation  whieh  was  used  in  the 
eomputation  shown  in  Figure  68. 


/  = 


(O 

In 


(D.4) 


Figure  68:  Line  Diagram  of  the  Process  to  Convert  Radians  Per  Second  to  Hertz 


3,  Moment  of  Inertia  Subroutines 

Modeling  the  extremal  dynamics  of  a  centrifuge  required  determining  its  mass 
and  mass  moment  inertia. 

As  shown  in  Figure  69,  the  Annular  Centrifuge  modeled  in  this  thesis  consists  of 
four  main  geometric  components.  A  hollow-cylinder,  or  dmm,  is  used  to  provide  the 
floor  the  centrifuge’s  occupant  would  perceive.  The  dmm  is  covered  on  each  side  by  a 
pair  of  end-caps,  or  discs,  that  serve  as  walls  for  the  astronaut.  The  centrifuge  is  then 
attached  to  its  prime  mover  by  a  shaft,  a  long  thin  rod,  which  transmits  the  rotational 
force  of  the  motor  to  the  centrifuge.  The  sum  of  the  masses  and  mass  moment  of  inertias 
of  each  shape  make  up  the  total  for  the  centrifuge. 
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Analytical  Term 

Meaning 

Units 

Jn 

Mass  Moment  of  Inertia  of  n-shape 

kg-m^ 

1 

Shaft  Length 

m 

m 

Mass 

kg 

r. 

inner 

Centrifuge  Inner  Radius 

m 

r 

outer 

Centrifuge  Outer  Radius 

m 

^shaft 

Shaft  Radius 

m 

t 

End-cap  (Disc)  Thickness 

m 

K 

Volume  of  n-shape 

m^ 

n 

Circular  Constant,  3.1415... 

Unitless 

P 

Material  Density 

kg/m^ 

Table  18:  Terms  Used  in  the  Mass  Moment  of  Inertia  Subroutines 


Figure  69:  Geometric  Decomposition  of  an  Annular  Centrifuge 
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a. 


Disc  Mass  Moment  of  Inertia 


This  VI  computes  the  mass  and  mass  moment 

of  inertia  of  a  disc  that  rotates  like  a  wheel. 

Radius,  m 

Mass,  kg 

^  0.00 

0.00 

Thickness,  m 

Mass  Moment  of  Inertia  kg  m® 

^  0.00 

0.00 

Mass  Density,  kg/m® 

y  0.00 

Figure  70:  Disc  Mass  Moment  of  Inertia  Front  Panel 

The  front  panel  of  the  Dise  Moment  Of  Inertia  VI,  shown  in  Figure  24,  is 
not  designed  for  direet  user  manipulation.  Rather,  this  VI  is  designed  to  be  a  subordinate 
VI  invoked  whenever  the  mass  and  mass  moment  of  inertia  of  a  disc  are  required.  The 
input  controls,  Radius,  Thickness,  and  Mass  Density,  are  linked  from  an  outside  source. 
The  necessary  computations  are  performed  with  the  results  passed  to  the  mass  and  mass 
moment  of  inertia  output  indicators. 

The  computation  used  to  determine  mass  of  a  disc  shown  in  Figure  71  can 
be  determined  by  multiplying  its  mass  density  (/?)  by  its  volume  (V)  as  shown  in  Equation 
(D.5).73 

=  (D.5) 

The  mass  moment  of  inertia  of  a  disc  rotating  like  a  wheel  can  then  be 
computed,  as  shown  in  Figure  71,  by  using  Equation  {p.€)P^ 

(D.6) 


^3  Beyer,  p.  129 

Meriam  and  Kraige,  Dynamies.  p.  630 
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This  VI  computes  the  Mass  Moment  of 
Inertia  of  a  Solid  Cylinder. 


b.  Hollow-Cylinder  Mass  Moment  of  Inertia 

The  front  panel  of  the  Hollow-Cylinder  Moment  Of  Inertia  VI,  shown  in 
Figure  72,  is  not  designed  for  direct  user  manipulation.  Rather,  this  VI  is  designed  to  be 
a  subordinate  VI  invoked  whenever  the  mass  and  mass  moment  of  inertia  of  a  hollow- 
cylinder  are  required.  The  input  controls.  Outer  Radius,  Inner  Radius,  Drum  Length,  and 
Mass  Density,  are  linked  from  an  outside  source.  The  necessary  computations  are 
performed  with  the  results  passed  to  the  mass  and  mass  moment  of  inertia  output 
indicators. 


This  VI  computes  the  mass  and  mass  moment  of 
inertia  of  a  hollow  cylinder  rotating  like  a  wheel 


Outer  Radius,  m 

r)  0.00 

Inner  Radius,  m 
0.00 

Drum  Length,  m 

jo.oo 


Mass,  kg 
0.00 

Mass  Moment  of  Inertia  kg*m^ 
0.00 


Mass  Density,  kg/m® 
^  0.00 


Figure  72:  Hollow-Cylinder  Mass  Moment  of  Inertia  Front  Panel 
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This  VI  computes  the  Mass  Moment  of 
Inertia  of  a  Hollow  Cylinder. 


Figure  73:  Line  Diagram  of  the  Process  to  Calculate  Hollow-Cylinder  Mass 

Moment  of  Inertia 


The  computation  used  to  determine  mass  of  a  hollow  cylinder  in  Figure  75 
can  be  determined  by  multiplying  p  by  Vrod  as  shown  in  Equation  (0.7)7^  In  this  case, 
two  volumes  must  be  computed.  The  first  volume  is  that  of  the  overall  cylinder  from 
which  is  subtracted  the  volume  of  the  hollow  part  of  the  cylinder.^® 

=  PVdru,n  =  P  )  (0-7) 


The  mass  moment  of  inertia  of  a  hollow  cylinder  (Jdrum)  is  found  in  the 
same  manner  as  its  volume.  Equation  (D.8)  is  a  modified  form  of  the  equation  for 
Cylinder, in  this  case  the  mass  moment  of  inertia  of  the  hollow  interior  has  been 
subtrated.78 

5  (D.8) 

c.  Slender-Rod  Mass  Mass  Moment  of  Inertia 
The  front  panel  of  the  Rod  Moment  Of  Inertia  VI,  shown  in  Eigure  74,  is 
not  designed  for  direct  user  manipulation.  Rather,  this  VI  is  designed  to  be  a  subordinate 
75  Beyer,  p.  129 


76  Meriam  and  Kraige,  Dynamies.  p.  593 

77  Meriam  and  Kraige,  Dynamics,  p.  630 

78  Meriam  and  Kraige,  Dynamics,  p.  593 
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VI  invoked  whenever  the  mass  and  mass  moment  of  inertia  of  a  rod  are  required.  The 
input  eontrols,  Radius,  Length,  and  Mass  Density,  are  linked  from  an  outside  souree.  The 
neeessary  eomputations  are  performed  with  the  results  passed  to  the  mass  and  mass 
moment  of  inertia  output  indicators. 


This  VI  computes  the  mass  and  mass  moment  of 
inertia  of  a  slender  rod  rotating  about  its  length 


Length,  m 
^0.00 
Radius,  m 
^0.00 

Mass  Density,  kg/m® 

^0.00 


Mass,  kg 
0.00 

Mass  Moment  of  Inertia,  kg*m® 
0.00 


Figure  74:  Slender-Rod  Mass  Moment  of  Inertia  Front  Panel 


This  VI  computes  the  Mass 
Moment  of  Inertia  of  a  Slender  Rod. 


Figure  75:  Line  Diagram  of  the  Process  to  Calculate  Slender-Rod  Mass  Moment  of 

Inertia 

The  computation  in  Figure  75  used  to  determine  mass  of  a  slender  rod  can 
be  determined  by  multiplying  phy  Vrod  as  shown  in  Equation  (D.9).^9 

mrod=pVrod=P^rH  (D.9) 

Beyer,  p.  129 
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The  mass  moment  of  inertia  of  a  slender  rod  {Jrod)  is  determined  using 
Equation  (D.IO)^*’. 

(D.IO) 

d.  Centrifuge  Total  Mass  Moment  of  Inertia 

The  front  panel  of  the  Centrifuge  Total  Moment  Of  Inertia  VI,  shown  in 
Figure  76,  is  not  designed  for  direet  user  manipulation.  Rather,  this  VI  is  designed  to  be 
a  subordinate  VI  invoked  whenever  the  total  mass  and  mass  moment  of  inertia  of  a 
eylindrieal  eentrifuge  is  required.  The  input  eontrol  eluster.  Centrifuge  Dimensions,  is 
linked  from  an  outside  source.  The  necessary  computations  are  performed  within 
subordinate  Vis  with  the  sum  of  their  results  computed  and  passed  to  the  mass  and  mass 
moment  of  inertia  output  indicators. 


This  VI  provides  the  dimensions  of  the  Centrifuge 
to  the  individual  Mass  and  Mass  Moment  of  Inertia 
Vis  and  calculates  the  sum  of  their  results. 

Centrifuge  Dimensions 


Total  Mass,  kg 
0.00 

Total  Mass  Moment  of  Inertia,  kg*m^ 
0.00 


Figure  76:  Centrifuge  Total  Mass  Moment  of  Inertia  Front  Panel 


Outer  Radius,  m 
^  0.00 

Inner  Radius,  m 

C)  0.00 

Shaft  Mass  Density,  kg/m? 

Cap  Thickness,  m 

^0.00 

t)o.oo 

Shaft  Length,  m 

Centrifuge  Mass  Density,  kg/m? 

'j  0.00 

yj’oroF”” 

Shaft  Radius,  m 

'j  0.00 

Drum  Length,  m 
(j  0.00 

80  Meriam  and  Kraige,  Dynamics,  p.  63 1 
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The  formulae  used  to  eompute  the  some  of  the  results  of  the  subordinate 
Vis,  as  shown  in  Figure  77  are  imbedded  into  the  wiring  diagram.  The  equations  upon 
whieh  those  eomputations  were  based  are:*i 


J. 


centrifuge 


components 


disc) +  drum 


rod 


(D.ll) 


By  extension,  the  mass  of  the  separate  eomponents  is  also  additive: 

^centrifuge  ~  ^  ^components  ~  '^^^disc)'^  ^drum  ^rod  (D- 12) 


Figure  77:  Line  Diagram  of  the  Process  To  Calculate  Total  Mass  Moment  of  Inertia 


Meriam  and  Kraige,  Dynamics,  p.  593 
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